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ABSTRACT 


III 


Among the various techniques suggested for low bit-rate 
coding, sub-band coding seems to be very promising. It fills 
the gap between waveform coding techniques and ‘vocoding 
techniques. It is also a simpler compromise to complex trans- 
form domain coders. 

In sub-band coding, speech band is partitioned into 
4 or 5 sub-bands and the sampling rate of each sub-band is 
decimated to the Nyquist sampling rate. Thus redundancy 
present in the speech is also reduced. The transmission bit- 
rate is kept low by implementing a coding technique such as 
APCM, which is most suited for encoding sub-band signals having 
poor sample- to-sample correlation. 

This thesis addresses itself to the study of sub-band 
coders. Various guidelines for selection of parameters of sub- 
band coders have boon discussed and suggosticais made, A sub- 
band coder for speech signals has boon simulated on the DEC 
system 1090 computer using FORTRAN language. Simulation studio, 
have been done for 9.6 Kb/s and 16 Kb/s transmission rates. 
Conventional SNR’s and segmental SNR’s have been measured I 

and tho dynamic range of tho coder has boon assessed by I 

varying the input level, 

ADM has been inpl omen tod as an alternative to APCM I 
coder to study whether simplicity and adaptability of ADM can 
be exploited for encoding sub-band signals, I 
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CHAPTER 1 


INTfDDUCTION 

Speech is one of the main aspects in which the human 
race can be distinguished from the rest. The ability of 
speech in human beings inspired them to develop better 
mea ns of c ommun ica ti on s , 

With rapid developments in the field of science and 
technology, there is always an ever increasing demand on 
existing communication facilities, Vtlith shift from analog 
to digital communication, in the second half of twentith 
century, the communication engineers have to ensure the two 
following fundamental aspects in all their designs: 

a) minimizing the number of bits which must be trans- 
mitted over the communication channel to convey the 
given informatic«i within given fidelity requirements 

b) ensuring that bits transmitted over the channel are 
received correctly in the presence of interference 
of various types and origin. 

The representation of the same amount- of information with 
fewer bits, has twofold advantages. Firstly, the load on 
the existing communication channels can be reduced and 
secondly, the data storage capacity can bo increased, at the 
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same time. This twofold objective can be achieved 

by low bit- rate coding. Low bit- rate coding requires exploi- 
tation of the fact that speech signals contain large amount 
of redundancy. This redundancy- removal or da ta-c empress ion 
could be used to decrease bandwidth, to increase rate of 
transmission, to control probability of error and to reduce 
average signal power, besides low bit-rate generation. 

In this thesis, a particular class of waveform coders 
called sub-band coder for lew bit-rate speech transmission 
has been reported, 

1.1 SPEECH DIGITIZATIOH ; 

Techniques of speech digitization are remarkably 
varied. Their difference depend mainly on the properties of 
Speech and hearing, which arc exploited in the design. These 
factors also influence the resulting bit-rate. 

There is always a need for various speech coding 
techniques ranging from high to low bit-rates for various 
applications. Of course, the specific need win depend upon 
bandwidth and quality requirements. The range of bit-rates 
for transmission might be seen from Fig, 1.1, Speech digi- 
tization techniques can also bo broadly classified from 
Fig. 1,1. 


I 
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The left, high bit- rate side in Fig* 1,1, represents 
coding techniques tliat try merely to describe the acoustic 
Waveform of the signal within given fidelity criterion. 
Those techniques do no further coding of the source. This 
particular class of coders is called Vfaveform Coders. This 
class of coders includes high bit-rato pulse code modu- 
lation (PCM), adaptive pulse code modulation (APCM), linear 
delta modulation (LDM), adaptive delta modulation (ADM), 
differential pulse code modulation (DPCM) etc. These coders 
avoid liie coder complexity and reproduce speech with a 
quality sufficient for commercial purpose. Important perfor- 
mance considerations in v/avoform coding are bit-rato, mean 
Squared error (MSE) of reconstructed signal with respect to 
the original,, dynamic rang^,, implementation complexity and 
ruggedness towards transmission errors. 

The right hand side of Fig, 1,1, represents the 
techniques which assume a speech generation model and trans- 
mit updates of the parameters of continuously changing model 
characteristics to achieve low bit-rates. These techniques 
include all analysis synthesis methods, the gamut of voco- 
ders, linear and adaptive predictive coding, format synthesis 
and computer synthesis from stored text. The coders following 
those techniques are called parametric coders. Though this 
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Figure 1,2: Speech Digitizaticai Techniques 
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class of coders provides low bit- rate but it is highly talker 
dependent, sensitive to background noise and expensive. 

This classification has been illustrated in Fig, 1,2. 
1,2 LOW BIT- RATE GENERATION: 

The information capacity required to transmit or store 
the digital representation is given by 

I = B.Fg = Bit-rate in bits per second (l.l) 

where, 

F_ - Sampling rate (i, e, sarapies/second) 

O 

B - Number of bits per sample. 

It is desired to maintain the bit rate as low as possible 
while maintaining a required level of quality. For a given 
speech bandwidth, the minimum sampling rate is fixed by 
sampling theorem, Thereforo, the only v/ay to reduce the 
bit-rate is to reduce number of bits per sample. 

The various techniques discussed in Section 1,1 try to 
reduce the bit-rates in their own ways. Adaptive delta 
modulation (ADM) coders provide satisfactory performance 
at bit-rates above 32 Kb/sec and adaptive differential pulse 
code modulation (AEPCM') coders are suitable for 24-32 Kb/sec, 
range [l]. However , increasing n^d is being felt 
for efficient waveform coders in the 7,2 to 24 Kb/sec. range. 
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This is to benefit from the possibility of digital speech 
transmission over voice band provided by switched telephone 
lines and mobile VHF radios [2]* 

Though, the existing analysis and synthesis techniques 
provide low bit- rate but these are very expensive. Trans- 
form coding (TC) schemes also achieve the low bit-rate by 
suitable linear transformation of blocks of speech samples 
with the decorrelation being achieved in the transform 
domain [3], But the transform coding schemes are also very 
complex. 

A particular class of waveform coders called sub-band 
coder' provides a compromise between sophisticated and 
expensive analysis syn •thesis techniques and high bit- rate 
waveform coding techniques such as AUA and adaptive differen- 
tial PCM (ADPCM) in 1,1 to 16 Kb/soc. range [3,4,5], Sub-band 
coder is also simpler compromise to complex transform domain 
coders, 

1.3 SUB- BAND CODING: 

Sub-band encoding uses a frequency domain analysis of 
input signal instead of time domain analysis. It is based on 
the partitioning of the speech band into sub-bands and 
encoding the sub-bands individually. This technique offers 
aftractive possibilities for coding speech economically at 
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bi-t>-rates in the range 7.2 to 16 Kb/soc. combined with new 
emerging hard ware technologies such as CCD's (charged 
coupled devices), sub-band coding offers a practical and 
economical means of obtaining good quality in digital speech 
coding at low bib-rates. The partitioning of speech band 
into sub-bands enables to control and reduce quantizing 
noise in the coding. Each sub-band is quantized v/ith an 
accuracy ( bit-aiiocation) based upon perceptual criteria. 

As a result the quality of coded signal is improved over 
that obtained from a single full band coding of the total 
spectrmii. 

In sub-band coding, the speech band is partitioned 
into typically four or five bands by bandpass filters. Each 
sub-band is then lowpass- translated to dc, sampled at its 
Nyquist rate, and then digitally encoded using adaptive PCM 
(APCM) encoding. The encoded sub-band signals are multiplexed 
and transmitted. At the receiving end, the data is demulti- 
plexed, decoded and the sampling rate is increased to the 
original sampling rate. The sub-band signals are again 
filtered by bandpass filters and then summed up to recon- 
struct the replica of the original signal. 

Credit for early work on sub-band coder goes to 
Crochiere, Other early contributions came from Flanagan 
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and Vi/ebber [4,5, 6,2l] ,Crochiexe, Fiangan and Webber 
simulated the sub-band ccxier on computer and also carried 
out hardware implementation. They also compared the quality 
of the sub-band coder against that of AEPCM and ADf'A coders 
[6], Goodman and Wilkinson suggested a modification in APCM 
algorithm, which reduced the transmission rate to 7,2 Kb/s 
[22], Crochiere and Sambur further brought down the trans- 
mission rate of sub-band coder to 4,8 Kb/s by using varia- 
ble band coding scheme [23], 

1,3,1 Comparison of Sub- band Coding v^ith other Waveform 
Coding Techniques: 

Sub-band coding technique offers attractive possibi- 
lities for coding speech econcsnicaliy at bit-rates in the 
range of 7,2 to 16 Kb/sec, Informal listening tests were 
made by Crochiere, Webber and Flanagan [5] to compare the 
quality of tiie sub-band coder with that of full band encoding. 

When 16 Kb/sec sub- band coder was compared with 
16 Kb/sec, AEPCM, listeners preferred sub-band coder in more 
than 94>J cases, 'When the bit-rate of AEPCM coder was 
increased to 24 Kb/sec (3 bits/sample), the listeners rated 
sub-band encoded sentence as having higher quality in 34>^ 
cases. This clearly indicates that the quality of the 
16 Kb/sec sub-band coder is clearly preferred over that of 
AEPCM at the same bit- rate. 
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When 9,6 Kb/sec- sub-band coder was compared with ACM 
coder at different bit-rates of 10,3, 12,9 and 17.2 Kb/sec., 
the listeners preferred sub-band coders in 96>^, 82;^ and 6l>< 
cases respectively. 

When bit- rates of AEPQvl and AIM coders were increased 
to improve their quality then 16 Kb/sec* sub-band coder was 
found to have a quality comparable to 26 Kb/sec, AEJPCM coder 
and 9, 6 Kb/sec, sub- band coder had quality comparable to 
19,2 Kb/sec, A eM coder. 

Therefore, for same subjective quality, the sub-band 
coder has about a lO Kb/sec, advantage over ADPGM and AEM 
coders, 

1,4 ORGANISATION OF THE THESIS: 

This thesis consists of six chapters. In Chapter 2, 
the sub-band coder-decoder has been discussed and the main 
simulation steps have been indicated. Chapter 3 deals with 
the various criterion for the selection of sub-bands and 
integer band sampling. In Chapter 4, various aspects of 
decimation and interpolation are presented. It also deals 
with the specifications of the lowpass filter used in the 
process of sampling rate conversion. Chapter 5 provides 
detailed discussions on APCM and Aim coders and decoders. In 



10 


Chapter 6, the simulation results have been given and 
conclusions are drawn on these results. It also includes 
suggestions for future work. 

Under appendix A, input files have been discussed 
and signal flow-chart for the simulation program has been 
shown. Appendix B presents signal flow-chart for the 
filter design. Appendix C presents listing of the program. 



CHAPTER 2 


SUB- BAND CODER - SYSTEM DESCRIPTION 

For digital transmission, a signal must be sampled 
and quantized. Quantization is a non-linear operation 
and produces distortion products that are typically broad 
in spectrum. Because of the characteristic of the speech 
spectrum, quantizing distortion is not equally detectable 
at all frequencies. Coding the signal in narrower sub-bands 
offers one possibility for controlling the distribution of 

quantizing noise across the signal spectrum. This also 

» 

realizes an improvement in the signal quality. 

In coding any waveform, there are two basic concerns: 
the sampling rate and the number of bits per sample. The 
product of these two quantities gives the bit- rate. To 
substantially reduce this bit-rate, without making any 
compromise in the quality, use must be made of the ccais iderable 
redundancy present in speech in the form of non-uniform 
amplitude distribution, large correlation between successive 
samples and the consequent non flat spectrum. 

The redundancy removal in sub-band coders is done in 
frequency domain by decimating the sampling rate of each 
sub-band to the minimum acceptable sampling rate i.e. 
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Nyquist sampling rate of the sub~band. The lowpass 
translated, Nyquist rated speech samples have very low 
sample-to-sample correlation# which minimizes the redundancy 
present in the speech samples. 

The bit>*rate in the sub-band coder is further 
reduced by using such encoding techniques which use the 
least number of bits per sample, 

2.1 SUB-BAND CODING: 

The simulation of the sub-band coder on DEC-10 computer 
involves the following steps: 

i) Selection of sub- bands based on equal contribution 
to articulaticai index (Al) and to enable integer 
band sampling, 

ii) Filtering the speech band into the sub-bands 
selected in stop (i). 

iii) Decimating the sampling rate of each sub-band to 
2f^ whore f^ is the bandwidth of ith sub-band. 

iv) Implementing suitable encoding technique on the 
decimated samples in each sub-band. APCM and AUA 
coding techniques have boon used. 

v) Multiplexing the data of the sub-bands and 
inserting synchronization data. 
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vi) Providing synchronization detector and demulti- 
plexing. 

vii) Decoding the data of each sub-band, 

viii) Interpolating the sampling rate of data in each 
sub- band to the original sampling rate, 

ix) Filtering the outputs of the interpolators with 

another set of bandpass filters identical to those 
used in step (ii) in the coder side, 

x) Summing of the filter outputs of all the sub- bands 
to give reconstructed replica S(n) of input speech 
samples S(n). 

Steps (v) and (vi) have been eliminated assuming that 
the encoded data is available as input to the decoder after 
undergoing multiplexing and demultiplexing with synchro- 
nization* 

The block diagram of the sub-band encoder is presented 
in Fig. 2.1, 

The functioning of the various sub-blocks of the 
sub-band coder and implementation of tho above mentioned 
steps have been described in the following paragraphs. 

An important element in the design of sub-band 
coder is the selecticaa of sub-bands that permit integer 
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band sampling and contribute equally to articulation index 
(Al). The speech band is partitioned into N sub~bands so 
selected by bandpass filters BPj^ to BP^j. The various aspects 
of selection of the sub- bands and integer band sampling 
are discussed in Chapter 3. 

The output of each filter in the coder is decimated 
to a sampling rate 2f^ where is the width of the sub-band 
and i refers to the ith sub-band. The decimators are 
denoted by DEC^^ to DECj^ in the figure. The decimated 
sub-band signals are digitally encoded by waveform encoders 
ENCODjl to e^COE|^* The encoded samples of all the sub-bands 
are time-multiplexed for transmission over the digital 
channel. 

At the receiver, the digital signals are demulti- 
plexed and decoded by decoders, DECOD^^ to DECOE^. The samp- 
ling rates of the decoder outputs arc interpolated to the 
original sampling rate of S(n), the input speech signal by 
filling in with zero-valued samples. The interpolators are 
denoted by INT^ to INT^^. The sub-band signals are reconst- 
ructed by filtering the outputs of the decoders with another 
set of bandpass filters identical to BP^ to BPj^ used in the 
coder. The outputs of these filters are then summed up to 
give a reconstructed replica 'stn) of the original speech 
signal S(n). 
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2.1.1 Choice of YJaveform Coders: 

As discussed earlier, the sampling rate of sub-band 
signals is brought down to the lowest sampling rate, adequate 
to describe the information content v/ithin the band. The 
next step is to select such a encoding technique, which keeps 
up v;ith the aim of keeping the transmission bib- rate low. 

Among the various waveform coding techniques available, 
adaptive pulse code modulation (APGM) has been found most 
suitable to accomplish the encoding of Nyquist rated sub-band 
signals having poor sample-to-samplo correlation, 

2.1.2 Input Signal; 

The actual speech samples have been used as input* The 
original speech samples arc sampled at 6 KHz and quantized 
to 11 bit accuracy. The bandwidth is 2940 Hz. The sentence 
used ia ’Have you seen Bill*? 

The digital speech samples are from Hewlett-Packard's 
FFT analyser of EE Department, Indian Institute of Science, 
Bangalore. First 2048 samples of the speech material have 
boon used, which correspond to tho segment ’Have you s..,.*. 
Tho segment contains a voiced section (’a*) and an unvoiced 
section( ’ s’ ) , Tho waveform of speech segment is shown 
in Fig. 2.2. The sampling rate of these speech samples is 
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increased to the sampling rate of the bandpass filters 
used in the sub-band coder. 

2.2 ADVANTAGES OF SUB-BAND CODING; 

The advantages of encoding in sub-bands over full 
band coding are listed as under, 

1, By coding the speech signal in sub-bands, quantiza- 
tion noise can be contained in the sub— bands to 
prevent masking of one frequency range by quanti- 
zing noise in another frequency range. 

2, By using separate adaptation for each sub-band, the 
quantizer step size can be adjusted according to the 
energy level in each sub-band* Therefore, sub-bands 
with lower signal oiergy will have lower quantizer 
step size and contribute loss quantization noise, 

3, The bit- rate assigned to each sub-band can be opti- 
mized according to the perceptual importance of each 
individual band. In lower bands, where pitch and 
forpiant structure must be accurately preserved, a 
large number of bits per sample can be used for 
encoding, where as in upper bands where fricative, 
and noise-like sounds occur in speech, fewer bits/ 
sample can be used. 
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SELECTICN OF SUB-BANDS AND INTEGER BAND SAMPLING 


The selection of sub-bands involves a variety of 
considerations, such as the number of bands, bandwidth of the 
sub-bands and their location. The next step in processing the 
sub-bands is to perform lowpass and bandpass translation 
before coding, A variety of techniques exist for performing 
lowpass and bandpass translations. However, one approach 
is particularly attractive for hardware implementation since 
it eliminates the need of modulators. It is based on integer 
band sampling discussed in Section 3,2, The sub-bands have 
to be so selected as to contribute equally to the articulation 
index and at the same time permit integer band sampling. In 
this chapter, these issues have been discussed and sets of 
sub-bands for various sampling rates have been proposed, 

3,1 CRITERIA FOR SELECTING SUB-BAI'JDS: 

A good compromise in the number of bands necessary 
for sub-band coding is generally found to be four or five 
bands, V/hen less than four bands are used, bandwidths become 
too wide and do not allow for full utilization of the 
advantages of sub-band encoding. Designs with more than four 
or five bands tend to consume bandwidths in transition 
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bands of filter in addition to requiring more hardware for 
practical implementation. 

A useful preliminary guideline for choosing sub-bands, 
is to partition the speech band into sub-bands that represent 
approximate equal contribution to tho articulation index (Al) 
under noiseless conditions. In this way each sub-band 
contains a significant portion of important frequencies of 
the sub-band, 

3.1,1 Concept of Articulation Index (AI): 

Articulation index (AI) is defined as a weighed 
fraction representing for a given speech channel and noise 
condition, tho effective proportion of the normal speech 
signal, which is available to a listener for convoying speech 
in toll igibil ity , 

Articulation Index is computed from acoustical 
measurements or estimates at the oar of a listener, of tho 
speech spectrum and of effective masking spectrum of any 
noise which may be present. From the articulation index 
measurement tests a curve has been derived plotting arti- 
culation index (Al) versus frequency [7,8,9]. This curve 
is presented in Figure 3,1, This plot is used for selecting 
the sub-bands which contribute equally to tho articulation 
index. 



amplitude in dB — Articulation Index 
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Figure 3,1: Articulation Index Versus Frequency- 



Frequency in KHz - 


Figure 3,2: Partitioning of Speech Spectrum into 
Sub- bands 
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3.1.2 Location of Sub-bands; 

Lower sub-bands should have narrower bandwidths and 
bandwidths should become progressively wider with increasing 
frequency. At lower bit rates, small gaps are permitted 
between bands to ccaiserve bandwidth and bit rate, as shown 
in Fig, 3,2. If the gaps are too large, the noncontiguous 
bands produce a reverberant quality in the signal. However, 
some highly useful compromises can be achieved between 
transmission bit- rate and quality. 

3.2 INTEGER BAND SAMPLING: 

Sub-bands are lowpass translated before coding to 
facilitate sanpling rate reduction and to realize any benefit 
which might accrue from coding the lowpass signal. 

The lowpass translation can be performed in a variety 
of ways. In the present work a technique named integer band 
sampling has been used. This technique is better for hardware 
implementation and it also avoids the use of modulators as used 
in other conventional techniques for lowpass translation. 

3.2,1 Sampling Theorem for Bandpass Signals: 

Signals with bandpass spectra can also be represented 
by their sampled values. The spectrum of a bandpass signal 
is shown in Fig, 3,3, 



iinimum Sampling Frequency 
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X(f) 




Figure 3.4: Minimum Sampling Frequency for a Signal 
Occupying a Bandwidth 
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If a bandpass signal x(t) has a spectrum of bandwidth 
and upper frequency limit f^,,, then x(t) can be represented by 

X XU 

instantaneous values x(KT£,) if the sampling rate f_ is 
2f 

where m is the largest integer not exceeding 
(Higher sampling rates are not always usable unless they 
exceed ). If the sample values are represented by 
impulses, then x( t) can be exactly reproduced from its samples 
by an ideal bandpass filter H(f) with the response 



< 1^1 < 


elsewhere 


(3.1) 


The sampling rate for a bandpass signal depends on the ratio 
If » I- then the minimum sampling rate 

approaches 2 A sketch of versus fg/B^^ is shown in 

Fig, 3.4. The exact reconstruction occurs when 


f,= (3.2) 

where m is an integer satisfying the oquaticn 

)-!<>"< W^x (3-3) 

Integer band sampling technique is based on this theorem. 
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3,2.2 Integer Band Sampling Technique: 

Integer band sampling technique is illustrated 
in Figure 3,5. 

The signal sub~bands are chosen to have a lower 
cutr-off frequency of mf^ and an upper cut-off frequency of 
(nH-l)fj^ , where m is an integer satisfying equation (3,3), 
and f^ is the bandwidth of the nth band. Typical values of 
m from 1 to 3 are most useful. The integer band sampling 
imposes the constraint that the ratio of upper to lov/er band 
edges of sub-bands be (raj^+l)/m^ where m^ is an integer that 
may be different for different bands. 

For hardware consideration, it is required that the 
sampling rates of sub-bands be derivable from a common 
cloclc. Furthermore, for digital or CCD hardware implement- 
ation, it is desirable to relate these sampling rates to 
the sampling rate of bandpjass filters by ratios that are 
integers. Finally, the requirements for multiplixing digitally 
encoded sub-band signals dictate that the transmission bit- 
rates of each sub- band be a rational fraction of the total 
bit-rate so that the data can bo framed and synchronized. 

Also a small fraction of this total bit-rate must be 
reserved for synchronizing and framing information. 







Resampled Signal 



Reconstructed Signal 


Figure 3.5: Integejo-band Sampling Technique and a 
Frequency-domain Interpretation. 
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3.3 CHOICE OF SUB- BANDS FOR VARIOUS SAMPLING RATES; 

The multitude of constraints discussed in the previous 
section greatly resticts the choice of sub-bands. To 
assist in the selection of sub-bands, it is useful to construct 
tables such as Table 3,1, It is assumed in this table that 
the sampling rate of the bandpass filter is 9.6 KHz. Column 1 
indicates the integer decijnation ratios that relate sub-band 
sampling rates to 9.6 KHz. Column 2 gives bandwidth, f^, 
and column 3 gives 2 f^ sampling rates for possible sub-bands, 
Coluim 2 through 4 specify choices for band edges f^^ 

(m^ = 1,2,3,,,. ). Therefore, all choices for sub-bands are 
discernible from the table once the sampling rate for the . 
filter is chosen. It is to be ensured that the sub-bands so 
selected from such tables should contribute equally to the 
articulation index (AI) by mapping the sub-bands on the plot . 
presented in Figure 3.1. 

The final choice of sub-bands is further restricted by 
distribution of bits per sample across bands, the total 
transmission rate and the multiplexing requirements. 

The transmission rate of each sub-band must be a 
rational fraction of the total bit rate so that the sub-band 
data can be multiplexed into a repetitive framed sequence. 

The lowest commcan denominator of these rational fractions. 
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Table 3.1 

Choice of Sub-bands for Integer Band Sampling and 
9,6 KHz sampling Rate, 


Decimation f. 

Ra tio 

2f^ 

3f^ 


1 

4800 

9600 

14400 

19200 

2 

2400 

4800 

7200 

9600 

3 

1600 

3200 

3800 

6400 

4 

1200 

2400 

3600 

4800 

5 

960 

1920 

2880 

3840 

6 

800 

1600 

2400 

3200 

7 

686 

1371 

2057 

2743 

8 

600 

1200 

1800 

2400 

9 

533 

1067 

1600 

2133 

10 

430 

960 

1440 

1920 

11 

436 

873 

1X9 

1745 

12 

400 

800 

1200 

1600 

13 

369 

738 

1108 

1477 

14 

343 

686 

1029 

1371 

15 

320 

640 

960 

12X 

16 

XO 

600 

900 

1200 

17 

282 

565 

847 

1129 

18 

267 

533 

800 

1067 

19 

253 

505 

758 

1011 

20 

240 

480 

720 

960 

21 

229 

457 

686 

914 

22 

218 

436 

655 

873 

23 

209 

417 

626 

835 

24 

200 

400 

600 

800 

25 

192 

384 

576 

768 


29 


Table 3.2 


Suggested 

Band Coder 

Sub- bands for 9.6 and 

« 

7.2 Kb/sec. 

four 

Sub- band 

Sub-baind 

Sub- band 

Edges in Hz 

Set No. 

Nos , 

From 

To , 

I 

1 

229 

457 


2 

457 

914 


3 

1067 

1600 


4 

1920 

2880 

II' 

1 

240 

480 


2 

480 

960 


3 

1067 

1600 


4 

1920 

2880 

III 

1 

253 

505 


2 

£05 

1011 


3 

1067 

1600 


4 

1920 

2880 
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including the fraction of transmission rate reserved for 
synchronization, determines the smallest frame size. 

Table 3,2 presents a few choices of sub-bands that 
can be used for 9,6 and 7,2 Kb/s four band coder. The 
gaps appearing between the sub-bands give the coder a 
Slightly reverberant quality. 

Choice of sub-bands for 10,67 KHz sampling rate can 
be derived from Table 3.3, Table 3.4 presents a few 
possible sub-bands derived from Table 3.3 for 16 Kb/s five 
band coder, 

3.4 DESIC3^ AND IIV1PLEMEI\ITATI0N OF THE FILTERS: 

The speech band is partitioned into sub-bands by 
bandpass filters, FIR filters of the order of 175-200 taps 
have been used. If wider transition regions are allowed, 
lower . order filters can be used at the cost of an increased 
reverberant quality of the coder, 

FIR filter is preferred over HR filter duo to their 
following properties: 

i) exact linear phase 

ii) no stability problem as ©ncountered in HR filters 

iii) realization is efficient 

iv) availability of efficient iterative design 
methods. 
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Table 3.3 


Choice of 
10.67 KHz 

Sub-bands for Integer 
Sampling Rate, 

Band Sampling 

and 

beclmaiicir 

Ratio 

^ f 


3fi 


1 

5333 

10667 

16000 

21333 

2 

2667 

5333 

8000 

10667 

3 

1778 

3556 

5333 

7111 

4 

1333 

2667 

4000 

5333 

5 

1067 

2133 

3200 

4267 

6 

889 

1773 

2667 

3556 

7 

762 

1524 

2286 

3048 

1 ® 

9 

667 

1333 

2000 

2667 

593 

1185 

1778 

2370 

10 

533 

1067 

1600 

2133 

11 

485 

970 

1455 

1939 

12 

444 

889 

1333 

1778 

13 

410 

821 

1231 

1641 

14 

381 

762 

1143 

1524 

15 

356 

711 

2133 

1422 

16 

333 

667 

1000 

1333 

17 

314 

627 

941 

1255 

18 

296 

593 

809 

1185 

19 

281 

561 

842 

1123 

20 

267 

533 

800 

1067 

21 

254 

508 

762 

1016 

22 

242 

405 

727 

970 

23 

232 

464 

696 

928 

24 

222 

444 

667 

889 

25 

213 

427 

640 

853 

26 

205 

410 

615 

821 

27 

198 

395 

593 

790 

28 

190 

381 

571 

762 

29 

184 

368 

552 

736 

30 

178 

356 

533 

711 



Table 


Suggested Sub~ bands for 16 Kt 


Sub-band Sub-band 

Set No. No. 


I 1 

2 

3 

4 

5 

II 1 

2 

3 

4 


5 
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Five Band Coder 


Sub- band 

Edges in Hz 

From 

To 

178 

356 

296 

593 

533 

1067 

1067 

2133 

2133 

3200 

190 

381 

333 

667 

593 

1185 

1067 

2133 

2133 

3200 
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3,4,1 Computer Aided Design of Linear Phase FIR Filter: 
FIR filters are characterized by their finite 


transfer function is a polynomial in 


> 


Their 


duration impulse response, b 


2 '' 


t • 


H(z) 


K 

E 

k=o 




(3,4) 


To design a filter of this type, one selects the coefficients, 
£hk^ , so that the transfer function has a frequency response 

4 A 

H(e‘^-'^), - Tc that approximates the design specifi- 

cations within certain tolerances. 


If the input signal is sampled at Fg, then the analog 
frequency f in Hz is related to the digital frequency 
through 

>v - 2itf/Fg (3.5) 

The desired frequency response H(z) may be specified cither 
in terms of real and imaginary parts, or equivalently in 
terms of its amplitude and phase, i.e, 

H(ej>^)=: E h. eJ>^ (3.6) 

k=o ^ 

It can be shown that, if the phase of H(e^^ ) is linear in >>,, 
the impulse response must be symmetric in the sense that 


*'k “ “V-k 


(3.7) 
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In this case, eqfuation (3,6) can be rev>/xitten as 

H(e^^) = [ E a cos(n>>)3 (3.8) 

n =:0 

for even K, where 
N - K/2 

ao a hj^ 

a 2 — 1>2,3,,,»,N (3,9) 

For odd K, equation (3,6) can be written as 

= [ S a cos (n-l)/2)>] e'J ^ 
n=o 

(3.10) 

where, 

N = (K+l)/2 ; a^ 2 hj^^^ (3.11) 

Leaving aside for the moment the linear phase term 
in equations (3.10) and (3,11), it is seen that the frequency 
response of the filter is given by a real cosine series, the 
coefficients of which are simply related to the impulse j 

respaase. The linear phase delay is determined by the length ! 
of the impulse response only. The problem of the design I 

of this type of filter beccanes, therefore, one of the 
finding the values hj^ so that the cosine series in eqn,(3,8)a: 
equation (3,10) matches the desired function of > as closely : 
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as possible. This approach, due to McCelian et.al. [lO] 
is very useful ‘ in computer aided design for a wide class 
of FIR filters, and the same has been used in this work for 
designing bandpass filters for the partitioning of speech 
band and as lowpass filters in sampling rate converters. 

The parameters of the bandpass filters are shown in 
Figure 3,6, The sub-band covers the frequency range from 
m^^ff to (m^+l) fj^. For practical reasons the filter passband 
must have a slightly narrower frequency range from m. f . + -^f 
to A transition region, Af# of the order of 

50 to 60 Hz was used in simulation with good results. A 
passband ripple of 0.173 dB and a filter stop-band attenuation 
of the order of 46 dB gave satisfactory results in simulation. 
Separate bandpass filter has been used for each sub- band, 
instead of bank of filters for the partitioning of speech 
band into sub-bands. This avoids the implementation 
complexity of the bank of filters, 

3,4,2 Implementation of FIR Filters: 

For implementing FIR bandpass filter, a non-recurrsivo 
realization has been used. In the non-recurrsivo reali- 
zation, the present filter output y^ is obtained explicitly 
in terms of only past and present inputs, i.e., previous 
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Frec{uency 


Figure 3,6; parameters of the Bandpass Filters, 



Figure 3,7; Filter Implementation. 
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outputs axe not used to generate the present output. The 
representation of non-recurrsive realization can be 
written as 


'n 


F( 


» 

n' 


m 

n-i 


(3.12) 


Fig. 3.7 presents the block diagram of filter realization. 
The discrete Fourier transforms (DFT) of input speech 
samples and filter impulse response were multiplied and the 
inverse discrete Fourier transform (IDFT) of the product 
gives the filter output. For taking DFT* s and IDFT s, 

FFT subroutinei of IMSL library have been used for 
faster execution. 



CHAPTER 4 


DECIMATION AND INTERPOLATION 

The process of interpolation and decimation aro 
fundamental operations in digital signal procossing-v The 
process of decreasing the sampling rate is known as 
decimation and increasing the sampling rate is referred 
to as interpolation in this context. 

In sub-band coders, the sampling rate of each sub- 
band speech signals is brought down to the Nyquist sub-band 
sampling rate for low bit-rate transmission. This operation 
is performed by the decimation process. For reconstructing 
the replica of the original speech from the low bit- rate 
rcpresv.nta tion , the sampling rate of the decoded sub-band 
signals is again increased to the original sampling rate of 
the input speech samples. This operation is performed by 
the process of interpolation. 

Two typos of computational issues generally arise 
in the design and implementation of decimation and inter- 
polaticn systems. The first issue involves the design of 
appropriate filter around which decimation or interpolation 
is based. The second issue involves the actual imple- 
mentation of the decimation or interpolation processing. 
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The design of decimation or interpolation filter 
involves the use of lowpass digital filters. Such filters 
can be designed in a variety of ways, e,g. window designs, 
equiripple designs etc. In this thesis, the computer aided 
design for designing optimum FIR linear phase filter [lO] 
has been used. This approach has been already mentioned in 
Section 3,4. 

The implementation of dec ima tors or interpolators 
involves the implementation of the digital filters in which 
the input and the output sampling rates are different. 

Because of this difference in sampling rates, a straight- 
forward implementation of FIR digital filter structure is not 
practical. Instead, special considerations must be taken 
in choosing the various specifications of digital filter 
and in efficiently implementing the digital filter, 

A general module of sampling rate converter has 
been discussed in the following section. The various speci- 
fications for lowpass filter to fulfil the requirements 
of the process of decimation and interpolation are discussed 
in Section 4,2, 

4.1 SAMPLING RATE CONVERTER: 

Fig, 4,1 shows the sajqpling rate converter, which can 
bo used to decimate as well as interpolate the sampling rates 
of sub-band signals. 
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Figure 4,1; Sampling Rate Converter 


Input Output 

Buffer Buffer 



Index Control of 
QBUF Update and 
COFS Selection 



Figure 4,2; Flowchart of Program for Sampling Rate 
Conversion. 
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The sampling rate converter converts the sampling 
rate of sub-band signals by factors of l/M where L and M 
are arbitrary positive integers. The input sampling rate 
f^ is increased by a factor L by inserting L-1 zero valued 
samples between each pair of input samples. Tlie insertion 
of zeroes is only conceptual and in practice, the input 
sequence 3 K(n) is treated appropriately, i.e,, it is assumed 
that behind each x(n) there are L-1 zero valued samples, 
when computing an output. This signal is then filtered 

with an FIR Ipwpass filter whose stopband cutoff frequency 
is f^/2 or L f^/(2M), whichever is smaller. The output 
signal of this filter is then reduced by a factor of M by 
keeping only one out of every M samples, i.e., only the Mth 
sample is computed. By the elimination of computing unnecessary 
output samples it has been shown in [l3], that the output 
y(n) can be computed by the relation 

Q*"l . , 

y(n) = S h(KL + (nra) .© l) x ( [ “ ] - K) (4.1) 

where <, > L implies the quantity in parenthesis modulo 

L and [ ] corresponds to the highest integer less than or 
equal to the number in brackets. The sequence h(n), 
n = 0,1,2..,.N-1 represents the coefficient of the FIR 
filter and is the number of taps in the filter such that 
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N < QL (4.2) 

and Q is an arbitrary positive integer. 

Fig, 4.2 shows the flowchart of the program which 
implements tho relation (4.1), Input data x(n) is supplied 
through BUFM and output data y(n) is received through BUFL. 
QBUF stores the necessary internal state variables and CX)FS 
stores the coefficients h(n), ICTR is a control memory 
which is generated by tho initialization program and is used 
to control the indexing of data and coefficients in the 
program. Tho program is based upon the one given in 
reference [15], The program performs decimation v/hen L is 
equal to 1 and interpolation when M is equal to 1. 

4.2 SPECIFICATIONS FOR LOViPASS FILTER: 

Tho purpose of decimating the sampling rate of sub~band 
signals is to obtain a reduced sample rate, adequate to 
describe the information content only within tho band of 
interest. To avoid aliasing at this reduced sampling rate, 
it is necessary to filter the original signal with a lowpass 
filter and then only sampling rate reductim could be 
achieved as described in Section 4,1. 

Similarly for regeneration of replica of original 
speech samples, the sampling rate of sub~band signals are 
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increased to the original sampling rate. Once the sampling 
rate is increased then it is necessary to remove the images 
of the signal spectrum that are centred at integer multiples 
of 2 tc/T, while leaving "the frequency below ^ unaltered. 

Therefore, whether it is a process of decimation or 
interpolation, lowpass filtering is required to avoid aliasing. 
Since it is impossible to realize an ideal lowpass filter, 
an ideal filter must be approximated, 

FIR filter is preferred in this case also as compared 
to HR filters due to the following reasons: 

1, An ideal interpola tor/decima tor has zero phase or 
almost linear phase. FIR filter can provide exactly 
linear phase. These filters are optimal in the 
sense that tho width of transition band between 
passband and stopband is minimum for given values 
of passband and stopband ripple and specified 
passband and stopband cutoff frequencies. Thus 
with FIR filters, the interpolation error due to 
phase non-linearity can be zero and error due to 
amplitude distortion can be made arbitrarily 
small, 

2, Though HR filters have recurrsive realization 
which serves computation time, the particular 


k 
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nature of sampling rate conversion problem makes 
FIR filters a better choice. 

In the following sub-sections, the specifications of the 
lowpass filters are discussed, 

4.2.1 Choice of Passband and Stopband Cutoff Frequencies: 

The choice of passband and stopband depends upon the 
Nyquist frequency (J2-) of the incoming sub-band signal. 

In passband i,e, 0 <. w ^ Wp, the frequency response 
should be close to 1.0. If sampling rate is increased by a 
factor of L, the filter gain must approximate L in passband, 
instead of 1.0. This is achieved by multiplying the impulse 
response samples of the filter by L [14], 

An error of + is allowed in the passband and the 
frequency response is required to be within + .v,g of zero in 
the stopband, as shown in Fig, 4.3. 

If the original sampling period is such that 
then (D- must be close to , as illustrated in Fig, 4.3, 
With the result, that the transition band < iji < is 
bound to be very narrov/. It implies that a large value for 
the length of the filter win be required* In such cases, 
it is reasonable to define only one stop band 

^ 


(4.3) 
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Fig. 4 



However in cases, where is significantly higher than the 

Nyquist frequency, the transition band between passband and 
stopband can be wider and it makes sense to define stopbands 
around §ach integer multiple of — ^ . 

4.2,2 Order of the Filter: 


The order of the lowpass filter mainly depends upon 
desired ripple specifications and cutoff frequencies. Given 
these specifications, the order of the filter can be 
calculated by the following relation [13,16] 



~ f( )^F + 1 


(4.4) 



46 


where. 


A F - Width of the transition band normalized to the 
sampling frequency. 

^co (Sp, %) - [5.309 X 

+ T.lKxio”^ (log^QSp) 

- 0.4761] - [2.66xi0“^(log^Q^p) 

+ 0.5941 (logio‘ap)+ 0.4278] 

(4.5) 


^ = 0.51244 logj^^C ^3/‘b55,)+ll. 01217 

(4.6) 

= tolerance in the magnitude response in 
the passband 

^ =s tolerance . in the magnitude response in the 
stop band. 


Generally, AF in (4,4) will be relatively small and the 
second two terms on the right side of (4,4) will be relatively 
insignificant compared to that of the first terms. Consi- 
dering these approximations, equation (4,4) can be written as 


D 


N *2 








In Table 4.1, the function , 

for some typical values of ^ ^ s 

P 


(4.7) 

Sg) is tabulated 
and for K=l,2,...4. 
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Table 4,1 

Tabulation of ^ 

Ss K=1 K=2 K=3 K=4 


0.100 

0.0100 

0.100 

0.0050 

0.100 

0,0010 

0.100 

0.0005 

0.100 

0,0001 

0.050 

0.0100 

0.050 

0.0050 

0.050 

0.0010 

0.050 

0.0005 

0.050 

0.0001 

0.010 

0,0100 

0.010 

0.0050 

0.010 

0.0010 

0.010 

0,0005 

0.010 

0.0001 

0.005 

0,0100 

0.005 

0.0050 

0.005 

0,0010 

0.005 

0.0005 

0.005 

0.0001 

0,001 

0.0100 

0.001 

0.0050 

0.001 

0.0010 

0.001 

0,0005 

0.001 

0.0001 


1.25 

1.46 

1.41 

1.63 

1,80 

2.02 

1.95 

2.19 

2.33 

2.53 

1.46 

1.67 

1.63 

1 .^ 

2.02 

2^25 

2,19 

2.42 

2.58 

2.82 

1.94 

2.15 

2.12 

2.33 

2.54 

2.76 

2.72 

2.94 

3.14 

3.37 

2.15 

2.35 

2.33 

2.54 

2.76 

2.98 

2,94 

3.16 

3.37 

3,60 

2.61 

2.81 

2.81 

3.01 

3.25 

3.46 

3.45 

3.66 

3.90 

4.12 


1.58 

1.67 

1.75 

1.84 

2.15 

2.25 

2.32 

2.42 

2.72 

2.82 

1.79 

1.88 

1.97 

2.06 

2.38 

2.47 

2,55 

2.65 

2.96 

3.06 

2.27 

2.35 

2.45 

2.54 

2,89 

2.98 

3.07 

3.16 

3.50 

3.60 

2.47 

2.55 

2,66 

2.74 

3,10 

3.19 

3.29 

3.38 

3.73 

3.83 

2.92 

3.00 

3.12 

3.20 

3.58 

3,67 

3.78 

CO 

ro 

4.24 

4.33 
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Referring to this table, the order of the filter is calculated 
directly for the tabulated values of ^ and "a . Equation 

Jr 

(4.4) has been realized in the subroutine for designing FIR 
filter in the simulation program. Giv^ the cutoff frequencies 
for stopband and passband, ^ and ^ , the subroutine 
calculates the length of the filter and designs the filter 
for the calculated length. 

The order of the filter, i.e,, N should be always 
odd. Because if N is even then linear phase FIR filter 
will have a delay of atleast ^ sample [l7]. This half 
sample delay itself corresponds to interpolation between 
samples, thus such a filter can not preserve the samples 
of the original sequence. So N should be odd. 



CHAPTER 5 


ViAVEFORM CODERS 

f 

Selection of a Suitable wavefoim coding technique that 
encodes and decodes the sub-band signals keeping the bit-rate 
low, and maintaining the required quality of the output 
speech is important. The various waveforms coding techniques 
available are PGM, APCM, logarithmic companded PCM, DPCM, 
ADPCM, LEM, AIM etc. 

Digital coding of Nyquist rated sampled signals is 
best acccMonplished using adaptive pulse code modulation (APCM) 
due to low sample- to- sample correlation of sub- band signals 
[ 1 , 3 , 4 , 5]. The adaptation logic and the various aspects of 
the selection of the coder, and the parameters of APCM coding 
for sub-bands have been discussed in Section 5,2, 

One of the objectives of this thesis is to study 
whether the simplicity and adaptability of ADM coders can be 
exploited for encoding sub-band signals as an alternative to 
APCM coding. The adaptation logic and various aspects of ADM 
coding have been discussed in Section 5.3. 

APCM coders and ADM coders fall under the category 
of adaptive quantizers and adaptive differential quantizers, 
respectively. The principles governing these quantizers have 
been discussed in the following section. 
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5.1 ADAPTIVE QUANTIZERS AND DIFFERENTIAL QUANTIZERS: 

In quantizing speech signals, on the one hand, it is 
desiied to keep the step size large enough to accomodate the 
maximum peak-to-peak range of the signal, on the other 
hand, it is desired to make the quantization step small 
so as to minimize the quantization noise, especially when 
the input signal amplitude is small. The amplitude of the 
speech signal can vary over a wide range depending on the 
speaker, the communication environment, and within a given 
utterance, from voiced to unvoiced segments. One approach 
to accomodate these amplitude fluctuations and maintaining 
a good 3^R throughout is to use a non-uniform quantizer. An 
alternate approach is to adapt the properties of the 
quantizer to the level of the input signal, in the adaptation 
technique, a higher a^R can be obtained as compared to the 
non-uniform quantizer. Adaptive quantizers can be classified 
according to whether they are slowly or rapidly adapting, i.o. 
syllabic or instantaneous. The APCM and ADM coding, used in 
this work follow instantaneous adaptation as proposed by 
Jayant [1,19], 

When adaptive quantization is used directly on samples 
of the input system, it is called adaptive pulse code modu- 
lation (APCM). The basic idea of adaptive quantization is 


• 4*1 
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to let tho stop size (or in gonoral the quantizer levels 
and ranges) vary so as to match tho variance of tho input 
signal. Adaptive quantization may be further classified as 
feed-forward and feedback adaptive quantizers as shown in 
Fig. 5.1. In feedforward adaptive quantizers, the amplitude 


Adaptive Cjuantizers 


Feedforward 

Adaptive 

Quantizers 


Feedback 

Adaptive 

Quantizers 


variable step 
size represent* 
tion 


Variable gain j 

representatiun Variable gain 

representation 


Variable step 
size 

representation 


Fig. 5,1: Classification of Adaptive Quantizer 

or variance of the input is estimated from the input itself. 
Whereas in feedback adaptive quantizers, the step size is 
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adapted on the basis of the output of the quantizer. The 
APCM coder used in this thesis follows feedback adaptation 
of the step size. This is illustrated in Figure 5,2. 

Differential quantization scheme is motivated from 
the fact that there exists considerable correlation between 
adjacent speech samples. This correlation is significant 
even between the samples that are several sampling intervals 
apart. It suggests that the signal does not change rapidly 
from sample to sample so the difference between adjacent 
samples have a lower variance than the variance of the signal 
itself. 

Figure 5,3 shows the general differential quanti- 
zation scheme. The input to the differential quantizer is 
given by the equation 

q(n) = x(n) - ■x(n) (5.1) 

which is the difference between the unquantized input sample 
and its estin^te. The estimated value is the output of a 
predictor system P, In differential quantizer, this difference 
signal is quantized rather than the input. 

The quantized difference signal can bo represented 
as 

d(n) = d(n) + e(n) 


(5.2) 




Figure 5.2: Feedback Adaptetion of Step Size. 




Figure 5,3; General Differential Quantization Scheme 
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where e(n} is the quantization error. The quantized diffore~ 
ncG signal is added to the predicted value xl^n) to produce 
a quantized version of the input, i.e* 

x(n} = x(n) + d(n) (5.3) 

Substituting equation (5.1) and equation (5,2) in equation 
(5,3) it is seen that 

x(n) = x(n) + e(n) (5,4) 

So it is quite evident that the quantized speech sample 
differs from the input only by the quantization error of 
the difference signal. 

The simplest application of the concept of differen- 
tial quantization is in delta modulation (DM). In DM system, 
the signal is over-sampled to increase adjacent sample 
correlation. At every sample, the sign of the difference 
between the input sample and the latest staircase approxi- 
mation to it is transmitted and the next staircase approxi- 
mation of input is increased or decreased by the step size 
depending upon the si go. The adaptive version of this 
technique is called adaptive delta modulation (ADM). 

5.2 AEAPTB/E PULSE CODE MODULAHON (APCM) : 




The choice of encoder parameter is determined in 
part by the static or laog term spectral characteristic 
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of the speech waveform. Fig. 5,4 illustrates typical long- 
term speech spectra (averaged over a sentence} based on 
measurements made in reference [7,18], Fig. 5.4 presents 
power spectral density of speech versus wraped frequency 
scale based on a constant (5 percent/division ) contri- 
bution to the articulation index (AI) [7], in order to 
illustrate the relative perceptual importance of the various 
frequencies. IWo possibilities for sub— band selection for 
low and high bit- rates have been illustrated. It is seen 
that across the entire speech spectrum there is characteri- 
stic drop in power density with increasing frequency. 

Across any one band, however, the drop in power density is re- 
latively small. Since Sub-band signals are low-pass 
translated and sampled at their Nyquist rate, they appear 
essentially as flat spectrum .signals at the low sub-band 
sampling rate and have essentially no sample- to-sample 
correlation. Figure 5.5 shows examples of sub-band sigials 
for band l to 4, Because of their low sample-to-sampl e 
correlation, encoding is best performed by adaptive PCM 
(APCM) [1,4]. 

5.2.1 Adaptation Logic and Coder Parameters for APCM Coder: 

APCM uses adaptive quantization directly on input 


samples. 
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Figure 5.4: Long-term Spectrum of Speech. 
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Figure 5,5; Typical V/aveforms of Uncoded Sub-band 
Signals 



57 


The step size adaptation strategy used in simulation 
for the APCM coders is based on one-word step size memory 
approach [1]. The adaptation logic is shown in Fig, 5*6, 


The coder input signal, denoted as for rth sample, 
is quantized to one of the 2^ levels according to the 
quantizer characteristic as shown in Fig. 5.7, where B is the 
number of bits per sample. The step size adaptation circuit 
examines the quantizer output bits for the (r-l)th sample 
and computes the quantizer size, o. ^ for the rth sample 
according to the relation 


r-1 




(5.5) 


Where, 


' ^ MIN ^ ^ MAX 


(5.6) 


and where is step size used for (r-l)th sample, 

a mul tipi i cat im factor whoso value depends on the quantizer 
magnitude level at time (r-l). It can take on 

one of 2^^ values, M^, , ,M2B~1, If the lower 

magnitude quantizer levels arc used at time (r-l), a 
value of M(Lj,_, 3^) = Mj_ loss than one is used to reduce the 
next step size. If upper magnitude levels are encountered, 
a value of M^^ greater than one is chosen. In this way, the 
coder continuously adapts its stop size in an attempt to 
track the short term variance of the input signal. The step 


is 
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Figure 5.6: Step size Adaptation Algorithm. 



Figure 5.7; Quantizer Characteristic for APCM Coder. 
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size is always kept between some minimum and maximum value 

'^MBI ^MkX* ^®®P®ctively. Typical values of for 

111 

4,3,2, 1^, 1^, ll' bit APCM coders are given in Table 5.1 

[ 1 , 63 . 

An interesting modification to the above algorithm has 
been proposed in reference [22]. This modification allows 
for encoding at an average bit rate of 1 +1/N bits/sanqjie, 
where N is an integer. In this approach the sign of the 
signal is encoded for each sample, r, and the magnitude 
of the signal is encoded with one bit every N samples. The 
step size adaptation is essentially that of equation (5,5) and 
the quantizer magnitude levels repeated for N-1 samples at 
the decoder. The sign bit transmits essentially the zero- 
crossing or phase inforimtion and the magnitude bit convoys 
the amplitude informa tic»i in the waveform at a reduced rate, 

\'iith this modification, sampling rate of the order 
of 7.2 Kb/soc* can bo achieved. 

5.2.2 Dynamic Range of APCM Coder: 

Tho dynamic range is the range of input variance, 
v/hich tho quantizer can handle. The quantities and 

in tho above algorithm represent practical constraints 
in the adaptation logic. Their ratio determines tho dynamic 
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TABLE 5.1 

APCM CODER PARAiliETERS 


B = 

4 

3 

2 li 

li 

^3 



Ml 

0*9 

0.85 

0.85 0.92 0.92 

0.92 


Mj 

0.9 

1.0 

1.9 1.4 

1 #4 

1.4 



0.9 

1,0 





^4 

0.9 

1.5 





% 

1.2 






M, 

1.6 







2.0 







2.4 






range 

that the 

coder 

can handle and 

their absolute values 

determine the 

centre 

of this dynamic 

range, ] 

For the 

input 

speech 

samples used 

in this work, a 

ratio of 

^MTIM 

= 100, 

was used, resulting in 

useful dynamic 

range of 

JViXiN 

about ; 

20 dB 

for the 

coder. 

The actual values of 


^MAX 

must 


be different for each sub- band, to match properly the dynamic 
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range characteristics of the sub-band codex to that of the 
long term speech spectrum. 

5.2.3 Allocation of Bits/Sample: 

A useful measure for assisting in the parceling of 
bits amc»\g sub-bands is the signal- to-quantizing noise ratio 
(S/N) as a function of frequency [6]. Fig. 5.8 shows typical 
S/N values as a function of frequency that are found to give 
preferred signal quality at bit - rates of 16 and 9.6 Kb/sec., 
respectively. At 16 Kb/sec, , it is found that good quality 
coding can be achieved with an allocation of 4 bits/sample 
( •: 19 dB S/N) in lower sub-bands, 3 Bits/sample (-,12.75 dD) 

in middle sub-bands and 2 bits/sample ('c 9 dB) in upper 
sub-bands. Similarly at the transmission rate of 9.6 Kb/sec 
allocation of 3 bits/ sample in the lowest sub-band and 

2 bits/sample in the remaining sub-bands is suggested, 

5.3 ADAPTIVE DELTA IvDDULATION (ADM): 

In order to exploit the simplicity of delta modu- 
lation at relatively low operating frequencies, adaptive 
delta modulator (ADM) has been proposed [l9,20]. In ADM 
coding the variable step size increases during a step 
segment of input and decreases while quantizing a slowly 
varying segment of input waveform and thus step size follows 



Sig:ial to quantizing noise 
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Figure 5.9: Adaptive Delta Moduiatic»i. 



63 


ti. lie- varying slope characteristic of the input signal, as 
illustrated in Fig, 5.9. 

5,3.1 Adaptation Logic for ADw coding: • 

nuj step size adaptation strategy used in simulation 
for ADIvi coder is based cn one-bit memory proposed by Jayant 
[ 193 . 

Fir, 5,10 shows the block diagram of adaptation logic. 
The AD/.i co<lor uses instantaneous , exponential adaptation in 



Fig. 5,10: Block diagram of AD.n logic 

the sense that the step size is changed at every sampling 
instant by one of two specific factors. The adaptation 
logic incorporates a one-bit meiuo3:y, that is, the inraediately 
past channel symbol is stored and is compared with the 







64 


incoming bit for a decision on the new step size 

Specifically, if the previous step size is denoted by 
the adaptation win be of the form 




m. 


= '’•Vi 

if 

Cr 

= '^x-1 


if 

Cx 

0 

H 

1 


(5.7) 


assuming that P and q are time invariant. 

Value of is given by 

Sgn '“r = S = V ' Wl) (5.8) 

Where, 

Xy is the amplitude of input signal and 

is the aiiiplitude of the latest staircase appro- 
ximation to this at the sa.iipling instant r. 


.. 


b,3,2 parameters of AOr. Coders : 

The crucial parameter of the adaptive delta modu- 
lator are time- invariant adaptation constants p and Q, 

The Smallest and largest allowable step sizes are other 
important parameters. 

In order to adapt to the signal during slope overload, 
it is necessary that 


P > 1 


(5.9) 
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and in order to converge to a constant input signal during 
a purely hunting situation (m = Q m , ) , it is necessary 
that 


Q < 1 (5.10) 

Jayant [19] has shown that P and Q should satisfy the 
reia tion 


PQ < 1 


(5.11) 


for stability, i«e., to maintain the step size at values 
appropriate for the level of the input signal. The values 
P = 1.5 and Q = 0,6566 have given good results in simulaticaa 
for wide range of input speech samples [l9]. 


Secondly, the step size limits should be chosen to 
provide the desired dynamic range for the input signal. The 
ratio of large enough to maintain a 

high SN’R over a desired range of input signal levels. The 
miniiiim step size should bo as small as is practical 

so as to minimize the idle channel noise. 

In the present v;ork, the ratio of has been 

used equal to 100. The values of adaptation constants 
P and Q used are 1.5 and 0.6666 respectively. 



CHAPTER 6 


SIJVULATim RESULTS CONCLUSION 

In this thesis, simulation studies have been carried 
out on real speech samples for 9,6 and 16 Kb/c cub-band coders. 
Four band and five band sub-band coders have been designed 
for 9.6 and 16 Kb/s transmission rate respectively. 

Combinations of four, three and two bit quantization have 
been used in these designs. 

The performance criterion used for simuiatim study, 
have been discussed in Section 6,1, Section 6,2 presents 
the results obtained from the simulation study. The details 
of the simulation program, input files and signal flowchart 
are presented in Appendix A* The conclusions have been 
drawn in Secticn 6.3. 

6.1 PERFORMANCE CRITERION: 

iwo perfornence criterion have been used in simulation 
study. One is the conventional Signal- to-Noise ratio (a^IR) 
and the other is segmental averaged asiR- ( SBGSNR) , 

6.1.1 Signal- to-Noise Ratio (34R): 

The long time averaged signal- to~noise ratio (SNR), 
in dB, is given by 
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mR = 10 logj^Q 


M 

E 

ifrl 

M 

I 



(Xi-x,)- 


( 6 . 1 ) 


Where i = l,2,..,,Mare the original signal samples 

and Xj_, i =s 1,2,..., Mare the corresponding reconstructed 
(processed) signal samples. M, the total number of samples 
processed, is 2040 in the present study. 

In measuring the input and output signals to the 
Sub- band coder, it is generally desirable to compensate 
for the effects of filtering in the coder [21]. This is 
done by an arrangement shown in Fig, 6,1, 



Fig. 6,1; Circuit for evaluating GNR of sub-band coder 
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The input speech signal 'S(m) is sub- band coded to form the 
output speech signal X(m), It is also filtered with the 
same filters used in tlie sub— band coder to generate a 
compensated reference signal X(m), which is used as the 
input signal in equation. (6.1). Thus, the 3Nia defined here 
is strictly a measure of coder distortions and is not 
affected by bandlimiting or group delay in the coder. 

6.1.2 Segmental SNR (3>IRSBG): 

This is another performance measure being used now a 
days. It is believed to be more representative of the 
subjective quality of the output speech signal than the 
average would indicate. Segmental averaged SNR, a^RSEG, 
is obtained by first calculating the SNR's (as defined 
in Section 6,1.1) for short segments and then averaging out 
these per-segment SNR's, 

Let i) be the SNR of the ith segment in dB, 

Then, 

1 L 

S'JRSBa = f' Z SNR(i) (6.2) 

^ i=l 

V#tere is the total number of segments processed* The 
segment duration is typically taken as 20 ms [2], A segment 
length of 128 samples has been used in the present simulation 
study. Then, L » » 16. To prevent any one segment 
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from dominating the average, the value of SNR(i) may be 
limited within a suitable range* The effect of threshold 
3'3R(i) and the percentage of total segments on segmental 
SNR has also been studied, 

6,2 SIMULATION EXPERIjVB^TS: 

The simulation experiments were conducted on 2 sets of 
sub-bands for 9.6 Kb/s and 16 Kb/s transmission rates using 
APCM coders. These sets are given as follows; 


SET No. 

BAND EDGES IN Hz 


From 

Jo 

1 

240 

400 

( For 9- 6 Kb/s transmi- 

. 480 

960 

ssion i9te) 




1067 

1600 


1920 

2880 

2 

170 

356 

(For 16 Kb/s transmi- 

296 

593 

ssion rate) 




533 

1067 


1067 

2133 


2133 

3200 
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Using these sots of sub-bands, overall signal-to-noise 
ratio (a^R) and segmental signal-to-nois e ratio (SEGSNR) 
were measured. These results arc disc'ussed in Sub-section 

1 • 

Thcj dynamic range of the sub-band coder was also 
measured by varying signal input level for the whole system 
and individual APCM coders. These results are discussed in 
Sub-section 6,2,2. 

ADM coders v/ere also used in the simulation, for 
encoding sub-band signals, to s-tudy, whether the simplicity 
and adaptability of ADM coders can be exploited for encoding 
sub-band signals. Tliese results are discussed in Sub-section 
6.2,3. 

6,2,1 Over-all and Segmental Signal- to-Noise Ratio: 

The two sots of sub-bands, mentioned in Section 6.2 
v;ere used for simulation studies. The sub-band set was 
Simula ted, overall and segmental .SInIR' s were measured. Signal- 
to-noiso ratio of individual APCM coder of each sub-band 
was also measured. These results are tabulated in 
Table 6.1 and Table 6.2 for 9.6 Kb/s and 16 Kb/s trans- 
mission rates, respectively. The segmental has been 
measured with 0 dB threshold of BVaii). Fig. 6.2 shows 
the effect of threshold of Si'IR(i) on segitientai SNR. 
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Figure 6.2s Effect of Thre$hold SNR(i) on SEGSNR 



Table 6,3 


Variation of Segmental SNR with £NR(i) Threshold 
for 9,6 Kb/s Four-Band Sub-band Coder. 


— 



Percentage o'f 
segments* used 

SEGSNR 

(dB) ^ 

-10 

93.75 

5.801747 

- 3 

87.5 

6.435395 

- 2 

81.25 

7.113296 

** 1 

75 

7.846380 

0 

68,75 

8.589637 

3 

62,5 

11.28653 

6 

56.25 

12.18344 

7 

50 

14.57595 

10 

43.75 

15.87721 


*Totai number of segments = 16 



Table 6,4 


Variation of Segmental SNR with SNR(i) Threshold 
for 16 Kb/s Five-band Sub- band Coder. 


S'iRCi) 

(dB) 

percentage of 
segments* used 

SEGS^R 

(dB) 

-10 

93.75 

16.79234 

6 

93.75 

16.79234 

** '2 

93.75 

16.79234 

0 

93.75 

16.79234 

4* 2 

87.5 

17.90 

+ 4 

87.5 

17.90 

+ 6 

87.5 

17.90 

4 * 8 

81.25 

19.6467930 

+10 

81.25 

19.6467930 


*Total number of segments = 16 
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Ficjurc 6# 3s Dynamic Itanye of Sub-band Codox. 



Figure 6#4s Dynamic flange of Individual APCM Coder. 
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T*^bie. 6.3 ^Tict ^ Toblo 6,4 show the variation of 
segniental liNFi with the variation in SNR(i) threshold and 
the corresponding percentage of segments used. 

6.2.2 Dynamic Range tosurements : 

The dynamic range of the sub-band coder and indivi- 
dual APC/.i coders was assessed by varying input levels from 
0 d3 to -30 dB and 0 dB to 4- 20 UB, The results are shown 
in Fig, 6,3 for sub-band coder and in Fig. 6.4 for indivi- 
dual 4-bit, 3 bit and 2 bit APCM coders. 

6.2.3 Performance of ADM Coder; 


As an alternative to APCM coder, ADM coder was used 
for encoding sub-band signals. The sub-band Nyquist rate 
(fg =a 2f|) was increased by 5,10,15 and 20 times, the average 
a^jR, was measured in each case, 'Iho results obtained are 
listed in Table 6,&, 


Table 6,5 

Performance of ADI>1 Coder 


S, d o. 


1 

2 

3 


Number of times Sub-band SE'IR (dB) 

Nyquist sampling rate 

incrAas^ed 


20 

15 


10 


9. 6398220 
6.7524877 
4.019430 
0.7624895 


4 


5 
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6.3 CONCLUSIONS: 

codirsg offers e no\i oitemative in speech 
cncoclijnt^ ^^hich fills the beiv/een wsvsforni cociinci 
techniques and voCoding techniques. More generally, it 
is rolatetl lo a class of coders knovJn as transform coders 
[3], 3ut i t is a simpler alternative to the complex trans-, 
form coders, it offers a quality that is significantly 
better than con vent icarjai waveforras coders at low bit- rates. 

The design of sub-band coder involves the consideration 
of a lar^ju number of paramo^ors and trade-offs. For many 
of these meters, no analytical means exist for choosing 
them in an optimal v/ay. In this thee is it is attempted to 
provide some guidelines and insight for selecting parameters 
of sub-band coders. These guidelines are based upon the 
computer simulation work carried out, 

Tht:* iollov/ing conclusicns can be drawn from the 
simulation study undertaken. 

a) Linear phase Fin filter is a better choice for bandpass 
filters, used for partitioning of the speech band and 
for avoiding aliasing in sampling rate conversion, 

b) Filter length of the order of 175-200 taps is sufficient 
for partitioning of each sub-band of the speech band. 
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c) Individual bandpass filters must be used for each 
sub-band, in place of a filter bank, for easier 
implements ti on , 

d) A transition width of the order of 50 Hz in bandpass 
filter is sufficient, 

e) A passband ripple of 0,173 dB and stopband attenuation 
of 46 dB are good compromises for filter tolerance 
specification, 

f) The Nyquist rated sub-band signals have very poor 

sanple-to-saraple correlation, hence their encoding 
is best accomplished by APCM coder. For the same 
reason, differential coders do not perform well in 
a sub-band coder system, 

g) Signal- to-n Oise ratio increases with number of bits 
allocated per sample, 

h) When maximum-to-minimum step-size ratio of APCM 
coder is ICXD, and the minimum step size is chosen to 
give maximum SNR, the quality of the coder remains 
relatively constant over a range of input levels of 
about 20 dB, which is about 10 dB loss than the range 
reported by Crochiere [6,21], 
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i) The overall for 9*6 and 16 Kb/s coders have been 
found to be 10,36130 and 13.917365 dBs, which are quite 
close to the SNR's reported by Crochiere [6,21] for 

the same transmission rates. 

j) The average values of a^IR for 4,3 and 2 bit quantizers 
have been found to be 19, 12.75 and 9 dBs, which are 
quite close to the SNR's reported by Crochiere [6], 

k) Adaptive delta modulation requires a sampling rate on 
the order of 10 times or more of the Nyquist rate for 
good performance. So it seems that it is not suitable 
for encoding sub-band signals, keeping the bit- rate 
low [12], 

l) The simulation study has been carried out with only one 
set of speech samples consisting of 2048 samples. For 
a better critical assessment of the performance of the 
sub-band coder, this study should have been carried out 
with different sets of speech samples consisting of a 
larger number of samples. The average’ results obtained 
in such a case would have been a better indicator of 
the performance of the system, 

6.4 APPLICATION: 

There is an increasing negd of efficient 
waveform coders in the 7,2 to 16 Kb/s range to benefit from 
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the possibility of digital speech txansmission ovex voice 
band pxovided by switched telephone lines and mobile 
VHF xadios. Sub- band encoder pxomises to fulfil this 
need. 

Other potential applications of sub-band coders axe 
in the field of narrow band communication, voice storage 
application, voice coordination on digital data lines and 
for secure voice communications by digital encryption and 
transmission over conventional data lines, 

6,5 SUGGESTION FOR FUTURE WORK: 

The following investigations may be carried out as 
follow-up of the present simulation study, 

a) Subjective assessment of coder quality is 

important. As such extensive subjective tests 
could be carried out to gather data to 

supplement the SNR figures documented in this 
thesis. Different sets of speech samples at 
various sampling rates could be used for this 
purpose. 

b) The transmission rate of sub-band coder can be 
brought down to 7.2 Kb/ sec, by implementing modi- 
fication to APCM logic as suggested in 
Chapter 5 [22], 
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c) The txans miss ion bit rate of the sub- band 
coder may be further reduced to 4,8 Kb/s by 
implementing variable-band coding scheme for 
speech encoding suggested in reference [23], 

d) Hardware implementation using a general purpose 
signal processing chip like INTEL 2920 can be 
investigated. 
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APPENDIX A 


SUB- BAND CODER SIAULATIOM 

Non-real time simulation of the sub-band coder has 
been carried out on DEC SYSTEM 1090 at IIT Kanpur. 

The simulation program has been written in R3RTRAAJ-IV 
language. Some features of FORTRAN 10 language have also 
been used. PORTRAN language is preferred because it is 
supported by most computer system and it supports complex 
arithmatic. 

The program consists of the main program and several 
subroutines and functions. The main program calls the 
Subroutine to proper order. The program is designed to be 
interactive. The flowchart of the simulation program is 
presented to Fig, A.l. A full listing of the program is 
given in Appendix C. 

The various input files have been discussed in 
Section A.l, 

A.l INPUT FILES: 

The various input data required by the program in 
various subroutines is provided by input files. There are 
three different input files; for designing bandpass filters, 
for sampling rate conversion on the coder side and for 
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sampling rate conversion on the receiver side. These files 
have to be amended after considering various issues of 
design discussed in this thesis, before processing sub— bands. 
These files are shown in Appendix A.l, The details of 
various data supplied by these files are discussed as follows. 

Input File for Designing Bandpass Filter: 

As an example first line specifies that a 175 tap FIR 
filter is to be designed with 3 bands and 0 grid density. 

The second line specifies that the desired amplitude response 
is i in pass-band and 0 in stop band. The third line specifies 
that the errors in passband and stopband should be weighed 
by 1 and 10^ respectively. Line number 5 to 8 specify band 
edges for the required passbands and stopbands. These input 
frequencies are normalized to 0.5, In the first line the 
length of the filter is indicated as 0, if the lowpass filter 
is to be designed. In this case, the subroutine calculates 
the length of the filter by equation (4,4) of Chapter 4 and 
designs the filter accordingly. 

Input File for Sampling Rate Conversion on the Coder Side: 

This file can be divided into segments of 7 line 

ft 

each. Each segment provides input data for one stage 
decimation/ Interpolation, If the sampling rate conversion is 
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to be done by a factor of more than 15, then two stage 
interpoia tion/decima tim is used. 

In each segment, the first line indicates sub-band 
number, the second line indicates values of M and L. M is 
the integer by which sampling rate has to be decimated and 
L is the integer by which the sampling rate has to be 
increased. Line number 3 to 5 provide the lowpass filter 
specifications as explained earlier. The last line of the 
segment gives the ripple in passband and attenuation in 
stopband in dBs, 

Input File for Sampling Rate Conversion on the Receiver Side: 

This file is divided into segment of 2 lines each. Th^ 
first line indicates the sub-band number and the second , line 
indicates the values of M and L, 
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APPENDIX A.l 

INPUT FILES FOR SIMJLATICN PROGRAM 


Input Files for Designing Bandpass Filters 


Line No « Input Data 


1 175 3 0 

2 0 10 

3 10 1 10 

4 

5 0.0 0.025 0.0302683 0.04479 0.05 0.5 

6 0.0 0.05 0.00552083 0.0947916 0.1 0.5 

7 0.0 0.1111438 0.1166667 0.16144583 0.16666 0.5 

8 0.0 0.2 0.2052083 0,2947917 0.3003125 0.5 


Input File for Sampling Rata Conversion on the Coder side. 
Segment No . Line No . Input Data 


2 5 1 

3 0 2 0 

4 10 

5 i 10 . 

6 0.0 0.0927083 0.1 0.5 

7 0.173 46 
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Segment No . 
2 


Line No . 

1 

2 

3 

4 

5 

6 
7 

1 

2 

3 

4 

5 

6 
7 

1 

2 

3 

4 

5 

6 
7 


Input Data 
1 

4 1 

0 2 0 
1 0 
1 10 

0.0 0.0177083 0.025 0.5 

0.173 46 

2 

10 1 
0 2 0 
1 0 
1 10 

0.0 0.0427083 0.05 0.5 

0.173 46 

3 

9 1 

0 2 0 
1 0 
1 10 

0.0 0.0482292 0.0555208 0.5 

0.173 46 
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Segment No « Line N o* 

5 1 

2 

3 

4 

5 

6 
7 


Input Data 

4 

5 L 

0 2 0 
1 0 
1 10 

0,0 0,092708 0,1 0.5 

0,173 46 


Input File for Sanpling> Rate Conversion on the Receiver side: 
Segment No , Line No , Input Data 


1 

2 

3 

4 

5 


1 

2 

1 

2 

1 

2 

1 

2 

1 

2 


1 

1 4 

1 

1 5 

2 

1 10 
1 

1 9 

4 

1 5 
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APPENDIX C 


PROGRAM LISTING 


The following pages contain a listing of the 
Sub-band coder simulation program in FORTRAN^ 



♦ ♦ 



^ MATH PRUGHAm for S»B sand COnKR, 

C THTS is the main program for StMni,ATTNG SUR-RA^U cnUFii. 
C the 'VARIOUS 5U8-R00XIHES ALONG HTi’H THEIR rNPUT/nuTpiJi 
C PARAN^;fFRS ARE EXPLAINED AS THEY OCCUR IW fuE PkOGRAM, 



CONMOR/A4/NFILT,HP 
COWMON/AS/HPl ,KK 
COMMCN/AR/NN, AVAL 

common/ao/shrur 

OI"EMjtuN »!Pf 5,1 ,HP 1 C5,2050),HP3C5r2050y, AV4L(S,7uSt>) 

REM* S’IRPH 

integer RP,HElT.r(5) ,TCnoE, answer, storecs) 

npc .fU'!iT=2I,4jEVlCE='USR'j 

90ij 

ROO 'u'^tNATOR TO NOTE JYOU ARE REQUIRED TO KFYIN VAHTunii 

1 ,AS A'3U when NUTIFIFU ON THE TERMINAL. TilTS IS DijNfc tiY 

1 ‘EuYINGiN TUOvSE DATA FOUnH WITilTI " ",45 REOUTKED.') 

'ioi 

SOI FuR'iATfAX , 'please KEYIM NO. OF SOB-OARnS DESTKFu IN TJiF OuPER') 

AC''LPr*,K 

c 

r K denotes "JUMHER UF sub-bands TN the suh-ija«d coder, 

c 

^02 

90? FUR lA rC4X, 'ENCOUlNG TECHNIQUF ?M 
T^DE ‘iOJ 

90l FuR 1 ATC 4 X, 'KEYIN *1“ IF APCM CnuTNG DESIRED HR KFYTm "0" 

I TJ- ADM coding DESIRED') 

ACCEPft,TCDDF 

no I* Dpal ,K 

TdP=«P . , 

CALo ETuDESdoP) 

Ti^u VD 1 

f04 Fjn !i\TdX, 'DO YOU WANT TO DISPLAY FILTER FREOUF.^C Y PLOT? ' J 

TY»E yUb 

fOS FijR lATf d, ' kFYIH "1” IF PLOT DESIRED , OTHERWISE KEYTN hD''*) 

Answer 

TFIhNSWeR.EO.oIGO to lOU 
call gRaRH(IRP) 

TYiRa. yOA 

906 FOD,4Ai (4X,'fReOUENCY PLOT IS STORED IN FILE VkOPLT . OAT ' J 

loo CALL FLriHPdBP,!) 

call SRAtK(xnp,u) 

SinHF(IdPi-t laKK 
DO iU laJ »KK 
HPY(I0P,I)«HP1CIBP»I) 

10 CtiNTTH'fE 

call CDOIHGC IBP, ICODE) 

TYPE 907 

fit roPMATCOXr'Dfl YOU ^ANT TO KNOW SU8-BAND 



cr o 


R 


,Fl?,6, ' Ho ' ) 


DESlRhn nxWc:RvjToF KiiYt' 

nCSu^l * f iSwFK 
PAfiU S'!k 

WKT fF(S ,<i(jQ)5[^jRU,j5 

f!’ij’'<>'^Tf 4X, 'atiB-tiArjn snr is :',Fi?.6,' ns') 

•H) lt> rat,Yw 

’»Ft lPP,11»4VALCSP,t) 

r‘fi’ifTf<:rK 

srat?:! ctttP) 

CAt.i. FMtiPClBPtO) 

Fi/isn('iarx«2i,Db:viCEa'nsK') 

Kia2(! iR;C.4f.b SMMtKt ,6^) 

T^Ph; #to 

FOPHA f f 4X, 'DO Yf)IJ WAMT TO KNOw UVERALu S'IR?*1 
TYPK ill 

FUFi*lAf f 4X, 'fCEYlB "t« IF SwR Dt;SIRh.n,ni‘HKR4TSF K^Yl 

ACCEPT*, AilSWEK 

Tl- (A»l3i4cP,Eu,U)Cn TU 300 

ny IP ‘^p«l,K 

TBP»Bp;AKa5rORE(IBP+U 

rui I^P TbI,KK 

PFt lTHP,n*HPi(IBP,I) 

roMTtRML 

CAI.P SPATEI(IBP) 

CaT.E r»/Tl -IPtlBPfO) 

CUN'i'f 

»(l«2«»|H/cAfii4 Sl)M(Kl,603 

CAliU S>IH 

fAfa. SRGS'JR 

CAti, pPISfIR 

CONflrtJJC 

TYPE 

roPNAf(4X, 'OPERATOR TO NOTF.:THF PROGRAM EXECOTlUM 
THE IHPUT AHO OUTPOT SPEFCH SAMPT-ES ARE STORED IN 
FIliES *'(|R5fl»DAT AMD FUR59,DAT RESPECTIVELY,'^ 

STOP 

END 


IS UVER 
FILES 


SOWROHTINE FlLDESl- 

THIS SyBROWtlME DESIGNS FIR BANDPASS AMD LUWPASS 
FIIiTER# for yariows sus^bands. 


C»R> 


nnn nnnn n n n o m m n ^ 




TiGHin-aHTa D^IUSITY pSKT tu 16 

P’l>nK(2*N'UM0S)«Bfl<(in FlUCii ARHA^^t.r'vvnu A'lO iJPi>'=:H hnGF':G 

Fnu each rahu 

FXCuRa "lUSJ-uEStHRD AMRLXTUnt HFaPuilBf'') AKPAY.Ti 1:=? 1 Tu 
PASGHh Jy A‘!u V IN STHFRANy, 

VirXCMijAyn j-*<TD F!I:jCTIU« ARRAY Tu p.ACH BA'JU • i'M 'Us 

'^JCTGHbfF nt 1 A*1D iO RESRT.TH HASSHAJ-iH ^NO ST'^FP 
nuTp'if: 

hRSPCJHSH A^JD FREglltNCY FtiFO'lbF Ht tul FIjuTGR. 


ripinni, /Al/UFS,-«IT,ALPHA,TEXT,UFCFG,«nkTy,PI 2 ,Ar,Pi;<Y,X,Y,<-Kr.) 

ruMfiny/A^/CUiU^T 

FuMfin.i /A V'lX # NOLP , ItlANH 

Cu'^Mnivl/AA/HFlLT,HP 

mMEYSlQH rKXTC202),AO(2n2) , ALPHA ».?U7),XC202) ,VX7G?J 
niPCM^TiPI '!C202),HPu(b,4P4),HFf5,l7, 4 04) 

OESC^232) ,GRTU(i2i2) , ^f0272) 

niMEHiJlUM FI)Ge( 20) rFXCtO) ,WTX(1 U) ,nLVlATC JO) 
niMENBTuM itHFGA(50),RESPAC60) 

IKTEnKP RP,FIL1 / FIL2 rCOUHT , nFIL I ( b I , NX , HOLP, TB A IIP 

nuHjJLC PRECESinN UHERA,SHMAK,i 5 nMAC,HFSPA,AT'iM,PiPPT.E,SVHTF, 

I FELTAF 

nu'UiT.e PRECISION PI2,PI 
nrjtiBLB PRECISION An,nev,x,Y 
PI2II6.2« SIRS 307 179586 
PI«J. 1415926535899793 

AKPAYb IFXT.AU, alpha, X,X,H APE DTHFiMSiOwFu NFMAX2F? 

THF Arrays UFS, grip AMO WT are UIHFNSinuFO 16(NFMAX/2+2J 

WFMAXaAOO 

PIT.l*20 

ir(BP.Fy*0)FHil»FlLl + l 
ir(OP.HE.O3COU’«T*0 

l»EAt;(FtLI #*) NFXLTCCOUNt) ,NBANDS,LGRTD 

tF (HFtLT(COHNT),GT.?IFMAX,nH,NriLX(CnutlT) ,Es2.3KALL FHPOR 
IFCimANOS.LC.O) MBANDSsl 


ORIO OKNSIfY IS ASSUMED TU BE lb UNLESS SPECIFIED OTHERWISE 


9 m 

to 

too 




'>h 


ir(LCfHlU-GE, 0 ) LGRIDaie 
Ja» 24 'NBAMOS 

RCADCFTlit#*) CFX(a>,J»l,NBANDS) 
REAOCriLt,*) CWTX{JJ,J*U 1 «»ANDS) 

IFCBP,E 0«03 so *0 100 
140 10 l»t#RP 
READ (FII,1#009)AAA 
rO^MWCAl) 
eoiWYi^'Oe ' 

CEOQKIO) 

tn‘i*'f'»WlCOUWK,M* 0 ) 00 to uo 

tiiSilEiv 
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' 

r ' I ' ' ’ 'V ' - ■. ! 
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tin 
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C 
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c 
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140 
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C 

C 

C 

c 


ISO 




OKAU(FTLl,*) F<TPPLE, AI’Tm 
r ' ( '*11 T l'«20 # HFV tCFs 'DSK ' ) 

P I !(»♦♦( Ri PHI, l-:/20) 

A r'*’.i = iO**(-ATT»'/2u) 

PT*'Pr.haCPlPPT,fc;-.n/CRTPPijP4.t J 
uni oc^iPpLK) 
nr^'fi*u!iunio(A TTN) 

^vrUPaRiPPLK 

P lF»F'f'L'atO,Ou530y»CRlPpLEf»23 + U.o7l1 -1*HTPPLF-U , 4 /6l I ♦ AiT 
4 . (-0 ,a026u*(RlPrLK*'t<2 J-0,F>041 ♦PirpCjh:«U,‘<!27'^ ) 

AT'*’i.ai1 *Oi2l7tO,S1244’l‘SVPXP-0,5l2 4'l»ATTrI 
nK'T.T AKst.f^GP C I ) -FOGF C 2 ) 

tlFTUT(ruiJUTJaRtFPLK/nLL’lAF-ATTN+(uFLTAF^*2)+1 
TFff^nur,4FU,TfCnU?IT) ,2).£0.OJWFlUTlCO[JN'*’l = wFlT/rfCnufU) + 1 

flKHJttsi 

i^'KnaO 

nlijni)aRrifa*(CUUWT)/2 
%IUf*t)awriT,T(CnUNT)-?1'NUDD 
HFnuRaMFTliTtCOnRT) /2 

Trfunun,Fw*l*AND,NFG,EO,0) ('t(FCNS=flFCF6+l 

SET fip peuse grid. the wumbek qf putwts Vi ghtu is 

CfTbTfcR GENGTH+l)»GRTD UEf»iSXTX/2 

CRTDtUaFuGEd) 

n£i.r*t4nRiu*NFCNS 

OEGE«0*5/0RGF 

.Iitl 

Gal 

LHAHOal 

puP«ronE(G+i) 

T£MF«GRinCJ) 

CALCULATE THE MAGNITUUP RESPONSE A*^fD THE HEIGHT 

F0?ICTinii OH the grid, 

RESIJ>*KFFCTRMP#rx,WTX,£BAND) 

WT(J)«WATe<TEHP,FX,WTX,LBAND) 

GRIDCj)»TEMP+OELr 
irCGRil>M),GT,FUP) GO TO ISO 
GO Xn 14S 
GRlD€0-l)«Fyp 

0EStJ-t}»EFrCriIP,rX,WXX,l46AND) 

*»TCa-l)*WATi<rfJP,rX,WXX,LBAWD3 

MAWOaLSAHDft 

IFCl»8ANC»,GT*iiAN0SJ GO TO ISO 
GRIOCJ)»«OOi<i»J ^ . , ■ 

GO f 0' 4 40' ■: '■ ^ ' 



'1(1 ) • iri *J ♦ALPHA (MFCUS) 

rni J?5> 

H(tlfrUojao,5fAlPHA(l)-0,25f alpha C 
'UHZ)aO.O « Ai.5i 

nu rn ^50 

H(t J«0,2'?*ALPMA(NFCNSJ 
nij ins 

'U‘I) ALPHA (NZ-J) -alpha ( hPC^^o+^-J ) ) 
H(riF:NS)aU.5*ALPHACl).0.?5*ALPHA(2) 
nu 350 lal,HFCNS 
■JpOCruU.jT + l ,I)aH(I) 

CO'JtTNUL 

|«0 

TF (MODU.EOtl.AHD.NEG.Eg.OlUal 
nu J 60 K«j#NrcNS-i 
T3fH»’CHS-K) 

MpOtCU'lHTtl ,KtNFCHS)aH(I) 
ruJ^TrHtrt; 

TFfMH.'JE.OlGO TU 370 

Hx*«ni,PtiriAfjD 

GQ TO IJiO 
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CUHTTafE 
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CONflHOE 


CALCUUAfR THE FKEQUEMCY RESPOHSE 

no 4t0 IKA «i»5y 

OHEGAdKA) » 0,OlOOyOOOO*{IKA-t) 

SUMAR » n.O 
SUMAC * 0*0 
no 420 OtK* 1,NFCHS 

S0MAR»SUMAR+H(NXK)»CDC0S(PI2»0MEGAClKA>4(fiIA-l)) t 
I nCO-S(PI2*UMB<3ACIKA)#CHriLt(COUNl‘)-MIK))J 

SUMAC » SUMAC + HCNTK) ♦(USINCPI2»nMEC;A (IKA) » (H XK-l ) ) + 
1 nSlH(Pl2*aMEaACIKA)»€NFlLfCCnUNT)-NIK))3 

Cij^f Igllg 

fF(«Oun*RU*U SUMARb8UMAR*H(NFCNSJ»UCQS(PI2*0MEGACIKA} 
t tCMrCMS«t)) 

irc»l000»ca*l.5 ^ SUMAC»SUMAC-KCMPCNS340SINCPl2*0rtEGAtrKA) 

RgSPACtKA3'' M^'D»tR.tCSUIlARf»2 ♦ SUMAC»^2) , ' ' 

' tfWACJKA)' f.v>OvQW«K?l''OC««S?AaTO , 





1 


<»y^ P’'JfHnTCl‘iA,F*>,3,qx,f9,3j ' 

Rf'.TUnf. : 

Fwn ; 

P’U^CTirif< ti:rF(TEMP,FX,1«TXrLRANu) j 

FU^JCTinN xn CALCUtiAXE TriE DEJ51REn ^<AC 4 tTf»DF REFpni^Sc i 

Ao A F'’[«CrTurJ OF FREOUErtCX, ! 

f 

r»lMhMi,TUM FXC51,WTX(5) I 

FFFaFXCl.RA»IO) : 

RETURiJ I 

FnH i 

FUMCTiOi^ ^ATE(TEMP,FX,WTX,T.BAimD) ' 

I 

FUHCTinu xn CAtiCnijATF THF weight FHWCXTuN as rt tHr^rt'TU:^ j 

OF FRE^HENCY 1 

nilENGfO'l FXC5),WXXC5J I 

>?A’^b»wTX(GOAWD) I 

refuph ; 

rijn I 

SUFK^uTIME error 

TX«'i I 

FORMAX( '*♦♦**♦*♦*♦♦♦ error IN INPUT DATA + ♦ ) 

STOP 
FxO 


SUFH'^UTINE RFMEZCEDGF^NBANDS) J- 

THTii JIIBRUHXINE IMPuFMFnTB THE RF'lFZ EXCHANGE ALGURITH*^ 
FOR ThF weighted CHERICHFV APPRUXIMATION UF a ruNTTuEUUfa 
function with a sum Of COSINES, INPUT TO THE SUBROUTINE 
ARF A DEMSF ORTD WHICH REPLACES THE FREyUENCY AXIS, THE 
DESIREO FUNCTION ON THIS GRID, THF WEIGHT FUNCTIOM ON THE 
CHID, the number OF COSINES AND AN INITIAL GUESS Of THE 
FXTERHAL FREOUESKIFS.THE PROGRAM MTNXHIZE’S THE CHEBICHFV 
ERROR M UETERNINIHG THE BEST LOCATION UF THE EXTREMAL 
FHEaUENCIES (POINTS OF MAXIMUM ERROR) AND THEN CALCULATES 
the COEFFICIENTS OF THE BEST APPROXIMATION, 





ftUBHOUTINe REMEZ(EDGE,NbANDS) 

common n%/ DES,WT,ALPHA,IEXT,NrCNS,NGRID,P12,A0,i>Eif,X,I,GRir) 
DIMENSION EDGE(20) 

DIMEWSIUN IfiXTC202),AD(202),ALPHAC?02)rX(202) >YC2a2) 

DIMENSION D1S(3232),<JPI0C3232),NTC3232) 

DIMENSION AC2©2)#PC20n,0C201) 

OOnSliE PRECISION PI2,DNUM,DDEN,0TBMP,A,P,a 
OOtlWE PRECISION ^ AO, W/XrY ^ 

ftte MAXInaM',, NOMftfeR,;Or IfEPAflOIIS ''OF as 







tJ 


! . ^ 


DhVba-t .0. 

«ZsuFcns+i 

Tf7,7Kr)rrrjs+? 

f'l‘’’i;pso 

100 ruHTTi^nt 

T t XT ( li 7 Z ) =H ijP 1 D 1 1 

f'flTt.PaMlTliR + t 

Ti- ffni'P'H.UT.ITPMAX) GU TO 400 
nij ItO jai.ME 
nTFMPaGKXDClKXK.T)) 

OTFHPsOCnSCDTfc^Mp^pTa) 

110 yi.T)»uTt;,Mp 

JbTs(NFCMS»l)/lt3 + l 
t»o no j*i,N7 

Jatal 

120 AU(0)aD(JJJ,MZ, JFI) 

nfj(fj.iao,o 

ODF«l«0.0 

Fot 

nu i3o j*i,HZ 
taTKXKj) 

OTF:MP*AOCj)»f)tSCL) 

DU»tH«OOUM+nTE(vlP 

f>TEMPaK*At>Cj)/WT(li) 

ddf:n*j)Oeh+otemp 

IJO Ka*K 

OEV*O0UM/0nEN 

TFtDEV.GT.O.O) NOa-i 

DEVa^KO^DEV 

Km^W 

no 140 j9lfHZ 
LalEXTCj) 

DTRf4P»K$neV/WTCL) 

Y(J)a0ESC!i)tDTEMP 
140 Ks-K 

rrCoEV.GK.OEVE) GO TO 150 
fAliO OUCH 
GO TO 400 
150 OEVtsoEV 

vlCHIiiGe*0 
Kl*ieXTtl) 
fCN2!*IEXfC?lZ) 

ICI|OW»0 
MaT**»NU 
J*t 

« 

SKARCH *^0R tHE EXTRE(4A0 FREOUENCIES OF THE REST 

APPROXIMATTON 

. lff4k'eE.,^fErW'-tO 300 ' ' ' .t'' , 

'KW»IEXfWtl' .. ' . 


210 




220 

225 




TFCJ.U0.2J li 
ritHpaorv 
TFfL.uF.fCUP) 
RnPanuPcr.^HZ) 


=CUMP 
Gn TU 220 


FHP*(ljRH-UF3CL))%WT(t,} 
D'fFf-l'^aNUTI'FHR-COMp 
TFfOTFMpaiF.O.O) GO TO 

rfjHPaiitlx^KRH 

r.ar. + i 


tFCL.GR.KUP) 

Ffm«Gfc;F(L,N2) 


GO TO 215 


FHR»CliRR*0ES(L3 JKWTCL) 
DlFMPaMuTfFKR-COMP 
tf(dtkmp.op*o,oj go to 

C0Mp«/jr?f*ERR 

no TO 210 
TEXT(J)»i:.-t / 


KG04»G-1 

.TCHNGE*JCHNGE'H 


GO in 700 

r.»o-i 


r.ari-f 


220 


215 


TF(i*,l,E,>CUO«) GO TO 250 
RHR«G6F(li,N2) 
ERR»(ERR-OES(L))t'WT(L) 
r)TEMR»flUf*RRR^CUMP 
irtOTEMp.GT.O.O) GO TO 230 

ircacHi^GE.Le.o) go to 225 

GO TO 200 

230 CUMP»«OT*ERR 
235 


TPfG,t,F.ICOO«) GO TO 240 
ERR«GEr(li,l«2) 
RgR«(ERrt»0ESCti))*WTCG3 
!>TEMR*WUT*«RR*C0MP 

ircaTE»«p*iie,o,oi go to 24o 

C0MP«WUT4ERR 
GO TO 235 

240 KliOWaifKfCO) 

IOT’Cd)«t+l 


250 

355 




JCRWGe«aCHRGETl 
GO TO 200 


C«i|ieiCTf4)tl 

ircaeHNe«i'GT*<^) gq to ^215 


fiStiTl > " , ; 

3waga?s‘s**, 

mmmsmm 


1 


760 

100 


110 


320 

325 

330 


340 

345 

150 

140 

no 

W' 'I ' > 


nu TO •'Uo 
Ta.T + l 

cr, in ?A>u 

Tff J.G'^'.nZZj GO TU J2U 
TF(M .r,T,XP:XT(U) KlslRXTtl ) 
TFfn.»lz.GT.TKXT(ij7J) KNZsTtXKNZ) 

HuTlBijllx 


I,»0 

KljOaKl 

roHP«sfMz*ct.ouooi ) 

r.uCKal 

TisTi + l 

tfcg.gj=:.kuoj on tu 315 

rKRaGEF:(G,WZ'l 

FRO»CEHH-0F.6(G) )*rtTCG) 

nTFi4PatJUT*KHR-C0Mp 

TF COTgMP.GPT.O.OJ GO TU 31U 

rij'’F«<^ot*eRK 

.|afIZ7 

nu in ?io 
riurK*6 
Gu rn 12R 

TK(GnCK.01,9) r,U TO 350 

ir(C0HO,GT.yi) yiaCo«p 

K1*1RXT(NZZ) 

t,aMORIO+t 

Kt,0rt*Kfl2 

ruMP*xt*ci,aouoi) 

trtU.GK .KLO^) GO TO 340 

ERR»GeRC0#M2) 

RKR»CERft*UESCL))*4»T{ti) 

OTR»PaUUT#€HR*COMI? 

TF(uTEMP,UR,0*0) go to 330 

JmHZ?. 


CaMP»liOT*ERR 

t.UC^*li0CK+10 

GO TO 235 

irCUOCK.RO.4) GU TO 370 
m 345 a»l#MPCWS 
t FX t ( WZ-J 5 » I K'XT C MZ-a > 


tgXTClIiiKl 

CO to too . 

KWlifEilTCMZZJ 
no 140 a«if»Ec»«5 
m%tunutti4*0 





1 Jh 


r 

c 

400 


405 


410 

415 

420 

425 

430 


$00 

505 

$10 

$50 



ru'J'A'TNifK 
nMi anFrfjfj-t 

riillal ,0fc.-06 

OTf-'HOa^HTUd) 

Xt*»/.7.i 3-2,0 
rws2*NFC’'fs«i 

f>KMf «i ,u/CN 
Ii*l 
F KK3(^ 

TrfKKK.EO.l) (in To 405 
r>TRHP»ncnSi:PT24GRXD(l )1 
r*NttMai}f'uSCPl2*GKlD(NGRID)) 

AA»2.U/(D'1TMP-ONOM) 

0B»-(aTi;;MP+DHUM)/CnTEMP-PNmO 

ruHTirtirc 

nu 410 jai,NFCM4i 

FT«(.T-l}*i)FLF 

XTa{)ruS{Pl2»FT> 

TF(KKK,fe;0,l ) GO TO 410 
XT»(XT-0R1/AA 
FT»AruS(Xl)/P12 
XE»XfGl 

TF(XT,G'f,XK) GO TO 420 
IF((XE-XTj,tT.FSH) GO TO 415 
G*G+1 
ao TO 410 

nu TO 425 

IF(tXT-XRl.Ll,FSH) GO TO 415 

aRTU(U«PT 

Ata)»GRECl,NZ) 

CUMTIHIIE 
TFCG.GT,!) GaL-l 
fUMTINUE 
GKiafUaSTEMP 
0I)FM»IPI2/CM 
00 5 m Jfcl'fNFCNS 
OTKWFwO.O 
Ol«l»M»Cvt-l)fBOgN 

GO to 505 
00 500 K*l#WMi 

OtRH*>»0fgllPTA€KTl)*OCOS(DHlIM4K) 

»TR«0«2 , 0*0TtMI»4X< 1 ) 

00 $»0 . ' ' < 

i»5| a 5:/CN 

li) /C:K 

C»f ^ ( NW 1 ) ' '' • ^ " 




1 ijl 


n(t IAC!,FCi^S-2)-ALr»hACMFrr,r>J 
DU bio 

Tf- f J*UT,;i,’l! ) GO Tu bib 
AAay.b*rtA 
PhaU,b»Hn 
515 Cg?»'lTaOfc; 

DC-Ut )3(J,U 
no b2u Kat , J 
n(» jaPCK) 

520 OCKia^.U'I'liD-fACK) 

r»(5'j*pf2)+A(1 )*2 i.0*AA 

roi b7b K»l,aMl 

525 Otlf)«p(K)+0(KHAA*ACKtl) 

JplaJ + l 

Do bio KaJoJPl 

530 DtK j|ai>(KHAA*ACK-n 

IF(J.K0.NI,'1t) GO TU 540 
ny bis K«i,U 
535 ')(K)*-ACK) 

0(l)«OCi>+AGPHA(NFCl^S-l-U) 

540 COHTTuDk: 

no 543 jal,MFCNS 
543 RUP<iA(J)»P(U) 

545 ronxTNire 

ir<Nrcws*aT,i3 keturn 
ALPHA (NfCWSf I 3*0,0 
ALPHAC5FCHSf2)»0,0 
MRffRC5,f 214) 

PfcTUPW 

Rwn 


nODBLE PRECISION FUNCTION 1 >(K,N,mj 


rOUCTlOM TO CALCULATE THE LAGRANGE IWTERPULATIflN 
COKfRlClEriTS FOR USE IN THE FUNCTION GEE, 


COWMOH /AI/ OeS#WTrAlPHA,iEXT,NFCHS,NGRID,Pl2|AD,DCif,XrIrGP£D 
DIWEMSIOH lEXf (202), ADC202 ), alpha ( 2023, X(202),y(202) 
niHEMSiON 0ESC3232),CRID(3232 ),wT(3232) 

DOUBLE PRECISION AO,DEV,X,Y 
DOUBLE PRECISION @ 

DOUBLE PRECISION P12 
DiiltO . ^ 

WCE) ' , ^ 

DO ,1 'W,M' ^ 

DO 'I 0»L^M>N ' ^ 


f'v 








’ " 






' ' ''*''C|***| ' 




1 <.p 


c 

c 

c 

c 


c. 

c 

c 

c 

c 

c 

c 

c 

c. 




t 

2 


FNn 


FRfXiS ION FIlhCTTuM f!eE(K,T;j 


TiGHll-o'' '“'’‘''■'f* HESPU'ISf HSIhG T,1- 

T.AUKAi, ,L IMTEKPULATinw FORMULA I»' THP '3 aP i:CL’'U‘P IC Fu*?h 


nr-inJ-Tn'J ,A(^0il,Y(.'>o?J 

PI LlaTOi f>KS(12l2),r,RTD(3232),wTt^232> 
nuPHf.F PRECISION P,C,0,XF 
nu'fui.i; precision eu 
nuRtifjlu PRECISION An,DEV-X,Y 
PaO.O ' 

XF»GRinCK) 

XF»PCUSCRI2*XF) 

nao,o 

ni) 1 

CaXF-ACJ> 

C*Auf3 VC 

DanfC 

PaPfC*YCJ) 

r,EF»p/p 

return 

END 

SUPKOUTINE ouch 
T lfRE I 

FORMATt FAILUHF TO CONVERGE V 

•OPHOliAOfaE CAUSE IS MACHINE PUUnDIWG ERROR'/ 

'OTHFR IMPtlOSE RESPONSE MAY «£ CORPECX'/ 

•OCHECK WITH A tREQURNCY HFSPONSFVJ 

RETURH 

END 




SURROUTIHE GRAPHCUPJj- 

THIS SU8H0UT1NE PLOTS THE FREOUEnCY- RESPONSE 
OF THE FILTER OBSIGHED M SUHROUTTNE FILuFSCaP^, 
INPUTl fRBOUENCY RESPONSE Of THE FILTER. 

OUTPIITl FREOUENCY PLOT. 




SOBROOTINE GRAPH(BP) 

OIMEWSIOW XC50),Y(50,I>,A(158I.IMAG4C5151) 
18TB0ER, BP^riLZ 
rtL 2 « 50 tBP^ " 

ORBW C OWjf :«f , OEVICEP/ OSK " , FILE* ' FRyPLT . DAT ' ) 
Spew t U'MI » DEVICE* '' OSK ' ) 

^'|»E«ICO«f* 2 «#DEVICB««''D$K%riLE*'FOR 25 .nAT') 







L,. 


! 


qyq 


C. 

c 

c 

c 

c 

c 

c 

c 

c. 


PuPMAi'r / 2 At) 

I i» UE V tCKx 'i)Sh * .Firii-'s' FkOl'T, f (jA'I’" I 

Ci.ni,KimitT.riw,nEncii='nsKii "'<'’‘'■‘•‘"'1 ^ 

Ct»'’> 5 E(Pi'< 11*25 , DEVICE*' nsK ' I 

RETKRri 

TADU UEHTiiT 

CADli 

E^n 


SUDKOuTiME FDTTHP(5p)j-. 

ImPuFMEnTJ the FTL'T’ER DESiGi^Fu .iY 
PURKnuTlHE FTL.nES(nt>) 

mlTPnit^TfliiSuLn^c^^^^^ filter IMPHeSE KEoPuP.'iF, 
nUTHJTJf TLTtRED SUB-BANU tilGMALS, 


SUDHnUTiHE FLTIMP(Bp,MnDFJ 
CUHM 0 N/A 2 /r 0 ’fMT 


cuMMnN/Aa/HXrNnGp^iaANn 
COMMOri /A42HFI1 :iT , HP 
roH,<tnfi/A5/Hpt ,ri 
EyniVAGEHCK(rMl) 

ttimnston XC2O50),HC2O5O),YC2O5o),TwK( 1 on), hPCP, 12 , 404 ) 

DIMENSIOH HPU5,2050) 

COMPUa CXC20S0),CHC2o5o) 

integer «P»MO0E,MfTLTC5),HX,NOLP,lRAMD 

DO 20 I»U*I 
X(i)»a,y 
H(t)»0,0 
CO!»«TiiOfWX*»fiP 
20 CO»«TIIIIJE 

DO 30 TwtfOriliKCOlUIT) 

30 COfITINfiB 

ir C WOE 0X50 TO 100 

DpBNCUMlT«47,OTVtCE«'DSK',FILE»'SPEl2Crf,DAT') 

Di) 40 tmipHih 

40 coMtiitie' 


100 


,4i; 

' 1.00 




00 'f 0 too 

iio«ao 4 ff 0 ri^i<T<oooi*t 3 






50 


t i 


CHCi )aCiir*UJtCHlI) ,0,) 

CO^'TlNUb 

CklA, FFT?CCX,MaWK) 
rkUL FFT2(CH,M,1WK) 
nu uo 1=1, N 
rx( i }j*cx(i)<-CKCi) 

60 rdriiifiiiK 

no /«» rsi,.\i 

XiaHL'AEi tCX( X ) J J X2»»AIMAG(CXt I ) ) 

CX( n»CMPLX(Xl,X2) 

70 CEJNA'lMilh' 

CALO frFH0R2CCX,M,IWK) 

TAlii. Kr'r2(CX,M,IWK) 

CAOb FFur>H2(CX,M,lWK) 

!»U HO Ia!,!\i 

y(n»Rt:AMcx(x)>/™AT(M) 

80 CUOIXMIF 

DO »0 

HP! (Rpa)»yci) 

90 COfllli^OF 

CtiOi>F(»riIX*47,FIL£«'’SPEi!;Ch.UAT') 

return 

END 

r 

iiiii aii i w ww <m iiiii»ii <ii i pi *w» iw i n i w m m i k" m - 

C — 

C SUBROUTINE SRATECBP, OPTION) 

C THIS SUBRIJUTINE IN CONJUCTION rtiTH SUBHOUTIDE SAoPCRP) 

C OECIBATES/INTEKPOLATES the sampling hates UF SUB-8ANU 

C ' SIGNALS. 

C 

C-. 

-'uwii^nniPi ■»»' I V - -mil .11 Hill .iiMii.iiuiiipwi,iipimwnn W i l l m il , | | I W > m $ i ip»i ( ii i iii | ii P ii i i i i>iiii iP W tw i<i>iwi i i ) iiii iii i i iii> p i i iiii i)Bii iiiiiii rifi i | i »ii i gB « i^ ^ 

SUBROUTINE SRATEtSP, OPTION) 

C0MM0N/A2/CQNNT 

CO«llON/A3/NX,NQLP,IfiAND 

COBNON/AU/rtiNT 


INTEGER CHANEL, BP, COUNT, NX, hOLP,IBAWD, OPTION, FIL3, HINT 
IF(0PTI0N.EQ,0) GO TO iOO 
riL3<»22 
GO TO 200 

too fll»l* 2 l 

200 CONtlNOe 

HiNf*0PtION' 

OPEN t , 0«VICE» 'OSS') 

REilNO rili3' — 

' eowwti^i ' - 

XBAN0#8P. ' ' . 


H ri>S"/ , ' «V, 


ircoR.c#.t'f.60;,to.ioo 

wW‘’ 



Mm: : 

®"'fp SOO' 




p 


roR'ai (/////, Al) 

nu i'D 6uO 

(FlhS,*) M,I, 
<JA«'101.Pi-IHAND 
CAUi. flLOdSCU) 

CAi.U t»A.4PtBl>,M,L) 
CltB.,TaCU!lfJT+l 

MUBP*P>ni^p+l 

Rt.AU tFir,'j,^,fc.NUs4U) CHANliT. 
IF IfHAtlEh.fiiQ.BP) GU TO 500 
Rfc:TURiJ 
EWD 


SUBKOUTlNt iiAHPCBF,M,L) I- 

TIIS antiRuUflHE DECIMATES/lMTfclPPObAXEiJ iHE PAMPl.T.^C 
RhTeS np SUB-BAMD SIGhAoS .IHXG SUBkDUTINj*, CAlLG SPl'UT 
TU irUTlAElZE AND TliEN CAGL& SPCOhV RUPPUnwG I/P OaTA 
TrtfnjlIGH BUFM AND TAKIi^G OUTPtll' DAiA FKOM tJ'IKL, 
INPUTlSUP-BA*^D SIG'NAD AND IMPULSE REPuNSE Dp THE UP FlLfEk, 
OUTPUTS DEC IMATEU/irifEKPULATED SUB-BANU SIGUAI.S. 


SUBHOUTIUB SAHP(BP,«,L) 

CUMMOR /SKCUM/ IQ,JQ,IL 
C0HM0r</A2/C0UNT 
COMMON/AI/NXrflOLP, IBAfJD 
COHHON/A4/NFILT,HP 
COMWON/AS/HPI ,NN 
eOMMON/Ali/HINT 

DIMENSIUN HPCS,12r404),hPl(5,2050l 

DIMENSION COErC200),COFSC350),0Bnt (t>50) ,TCTH(40) 

DIMCNSroM BaPL(J00O),BUrM(30O0) 

INTEGEH 8P#COUNT,«PILT(5),NX,NnLP,IbANl),HlNI 

IJ>»U 

ja»o 

ILi*U 

IFCMINT,BQ.0)G0 to too 

DO 10 K«(CNFILT(C0UNT)-i)/2+l)rNFlLTCC0UNT) 
COEFCK-(»FILT(COU«T>»l)/2)aHPCCUllNT4*l ,NX,K) 
C0«Tti4UE 
GO to 200 

DO 20 K»CCNFIt.T(COUNT)-l)/2+l)f8FlLTCCOUNT> 

COEf CK-C.irttiT(COONT)^t ) /2)*HP(C0UNTtl , f«X,K) *0 
CONTI WOE 
CONTINUE; ' v 

iNlTlilOIW CONVERSION ROUTINE. 






t J ? 


♦rlCHR 

DU il> Isl^i^N 
30 CUDiIMJe; 

c 

C DHOCK&S DkfA 

C 

MDKtM 10.4-1 )/M 

CADb SHCnf|V(BUFM,bUf L,Nl),g8UF,CuFci,lC'i*RJ 

;b<=ND»b 

DU 41! tal,NN 

Hl'l lHRa)«BUFD(I) 

40 CUNllMit 

return 

FND 



SURHnoTCNEiSRINXT 

TWlTIALIZAflUN t"OH SRCORV WHICH CUNi^EKTo tHe. dAwPDTDG 
RATE tJF \ SIGtIAL 6 t THE RATIO OK L/M 


•SUHOrjUTifie SRINlT(M,L,QBUK,N 0 fCDEF,N,COFS,NC,lC'i'R,Hi,lEHR ) 
COHMOH /SHCOM/ lQ,Jy,IL 

DIMENSIUN 9 BUF ( I ) ,COEFC 1 3 ,COFS( 1 ) , TC Trt C 1 ) 

M*OfcClMATION KAflO 
E*Ii<TEHKr)bAllON RATIO 
0 HUM 5 TATE VARIABLE BUFFER 

«g»SIZE or GBUF, GREATER OR EQUAL TO 2 »CTHE WEXT GREATEST 
INTEGER or M/L 3 

CUEr**ARHAX OF COEFFIClENTa FOR FIR INTEkPOLATING FILTER 
M«NO,UF TAPS in FIR INTERPOLATING FILTER 
COrSaSCHAMBLED COEFFICIENT VECTOR GENERATED BT SRiNlT 

nc«bi^e or cors vector, equal to or greater than 

L»<tHE NEXT GREATEST INTEGER OF N/L) 

ICTR»»CONTftOL ARRAI GENERATED BY SRINIT AND USED R 3 ! SPCORV 
NI*SIie OF ICTR VECTOR EQUAL OR GREATER THAN 2 *L 
IERK«*ERMOR CODE FOR DEBUGGING 

«U MO ERROR FOUND IN INITIALIZATION 
n aSUFCNQ) TOO small 
•2 cors CNC) TOO SMALL 
*3 ICTR CNI) too SMALL 

ieRR»D ’ ' ' ■ ' , ' ■' ' 

IfcWL ' ' ' ^ ^ 


COMPOfE I’g. 




IF tiiC.l.T.J'lF) Ifc.riK=2 
rJCF=(HFU/2 
F'ti s U 
C 

C UUT VBUF 

c 

DO il.‘ 

0 BHF(U s 0 ,y 

to r.ufifinJt: 

C 

C SCl'AMhfih COEFFICIENTS 

C 

I 'a I 

DU JO MU s 1,L 

DU 20 My a 1,IQ 

MX » (MU-l) + 

TF CMX.yT.NCF) MM a NCF - MX 
It (MX.GE.NCF) MM * MX-(N-NCF-U 
IF IMm.LC.NCF) COFS(I) a COEF(.MM)*FL 
TF (Mf4,UT,NCr) COfS(I) a 0. 

1 « lil 

20 CUNilNIlE 

30 COHTIfiOE 

C 

C SET OF MOVING AODHESS POINTER 

C 

jy » 10 

c 

C HENFRATE control ARkAT ICTR 

c 

EM » L*M 

ir (Nl.i,T,C2»li)3 IBRR * 3 
I.C «* 0 
MC * 0 

incr * a 

k a I I 

DO s»0 r » I, DM 
IF ILC.eO.O) XNCR * INCR+l 
IF <MC,liT.CM-i5) GO TO 40 

c 

C NO OF SAHFkES TO UPDATE QOUt 

C 

ICTRCKJ • IRCR 
IWCR • 0 
K • KFl, ' 

■■C ^JtARttit' MAftOH !»' cors VECTOR 

.€ ' ' ' ' , ’::,r ^ \ 


5U 


I 14 


MC * + 1 

CliNlTNlIt 
HtTuHN 
KNO 


C 

C 

c 

c 

c 

c 

c 

c 

c 



su^nnuTiMfc s srconv 


CUMtfP;KTi» 'IHE SAMPIjING RATE UF A SIGfJAL PY ThK RATiH L/'U 
SHItil'i MUST BE CALLED RHIUR Tlj CALLIML IHIS HUUiTRE, 


SUBhOUTlME SRC0NVCBUF«,bnFL,Nl),^PUF,CUFi3,xCtR ) 
common /SRCUM/ lQ,Ja,XL 

niMUHSlON BUFM(l)fBUFL(n,OBUt ( J ) ,COFSC 1) , ICThC 1 ) 

C 

C BUFO a INPUT DATA BUFFER UF SIZE nD+M 

C fttiFL a OUTPUT DATA BUFFER OF SIZE NO^'L 

C fU) a any POSITIVE INTEGER 

C QBUif a STATE VARIABLE BUFFER 

C COFS a SCRAMBLED COEFFICIENT VECTOR GENERATED BY SRINiT 

C ICTR * CONTROL ARRAY GEnEHATFU BY SHlNlT aNu .(SeD SRClUV 

C 

MB « i 

LB • 1 

L « XL 

^ DO SO 1 a 1,RD 
C 

C Me * INDEX ON flUFM 

C COMPOTE L OUTPUT SAMPLES 

C 


K • 1 - 

DO 40 0 » l,L 

JD • ICTRIKX 
IC » iCTR(Kfi) 
K ■ K+a 


c 

c 

c 

10 


'P::' 

f#; 


UPDATE OSOF 


IF CJO*K0.O) GO TO 20 
OdO#€JO) M iOFMCMB) 
agt •' JO #",10 
OftOrCJDl) » 00F«(CHB) 
Mi V • 

Jfi « J 0*1 

If cj 0 ,*M 4 or'J ,0 ».,-i 0 

'i..' ' 



t:’'t 


:t»^W ANR, store 




V ,1 ,1 r, 





I I s 


!H iO = l,i,j 

JC't* = + IC 

= Kt^( f Jtj 

r>u 1 = su^ + Q-iiiK( tQu)fCuF,s( j 

?0 Cii*! iltJOt. 

nuFI.tuHj = SUM 
Ui.i ” l,i'^ +• 1 

40 Cu'JT tilth, 

50 ruMTI-liJu 

tiLTuR.ii 

FriP 


C. 

c 

c 

c 

c 

c 

c. 


SuUKnirt'i'jf, C'<LluG;- 
TfCiS S j'jiih i!J I litift; 
'ihCiruOuh ni'Thf; 


CAhi^fj riiF sue. .'''.n £‘'11, t<u.< fht, 


FOR Ih ific 


th PhPU 


• H' 


r.'i'*-; III I t 4 f: t:nni.)(;tBp,Afisi'-t.pj 

PfAu VAU UlW ^TSiZt;, 

TMl l.r.hU B.>, A;)S«Uft,C'r.JU 

IF ) GO T(| 100 

ChI'G uI'hT{:chp,u 
CAGG VAhUh(tH'.J 
fAfOi aP IC’KOO) 

Tnl i.l 2 u<; 

loo CAhh VAiUANtOFJ 

Chhu hPCIC'KHP) 

fAhU A«’C -lOECHP) 

200 CiAJtM'lt 

Avh I’UHu 
f:ao 


c- 

c 

c 

c 

c 

c 

c 

c 

c 

c- 


r>j»^AnuTi.'jc, vahiamcbp) 

^•HiS S.JHOnUTlWli C?iLCyLATe.S FfeiA A , V AkIAUCf: A.jL' Si A.iMA'h' 
DRVCATIUiV UF TUK TIME SERIES OF SPEfiCR SAiiPOtS Ai' a-ht. 
£\piJT TO The encoder, 

IflPUTtOtClMATED SUB- 6 ANU SIGNAL, 
nuTPUT!SUl 5 -bAND SlGuAU 'MlHH UnTT VAkIANCF. 


SUPHOUTiitlE VAKlANtep) 

CUMi«lUfJ/A5/HPl,K 

CUMM0N/A6/STDEV 

dimension X(2050),HP1C5,2U50) 

REAL Sim,AMEAN,VAH,S 2 ,STDEV 

integer bp 



( if- 


'TPKuC j:iITsfiu,i)i;:vTCi:='i,n$*N' J 

'1p.’:.*(uru fsftl ^DF.VICC-'USK" J 
Sil 1 = << 

:>'J i<' r = l,K 

(,p, I ) 

XUJ-iHKpp, o 
3 J isn.j ! + Xtt) 

10 roaxii'i).; 

no -it' j=i,K 

S2=U(.r)-AMt^'\ O 
'/A'^ = '/APi’.'52 
20 CuOTT;,i 

vc ;-i ) 

r, rn<:vr:.'>.,<'ur vAc) 

ill.'! .!”•■( ol , f )h 
Mli .ii' I 

X(T. j=\f t)/.' i'OilV' 

r'H'’ U 'f ) 4( i) • 

10 Cl I 'i .‘Xu IK 

I’l/'o'U lu X t’=w'l r.CCs'OoK * X 
ni i.^FU ,nKi/iC-:='n;iK*) 

n xn 

C S.JiiHDuTlX!-: A0.1Ci)(0?») 

C I'ilis .'SliAK'lij f l.iK PK(XFuR;i,«i AoM CiilU'i'.i A''n iX-'C''*.; X .,C ,ic’ 

c sL'J.'iaks, 

C rNPlir;.>UR-'iAMU .'JlWOAi.. 

C TtjVt'iri;cnOF,i» A;'IU iJCCUUKU i;i!A-i:iA,4r) .'UF.lAi,. 

c 

StIf<K^Tl'Iii’'Aolcu73p7 

ncMi ;>A.-lPur-: , SS , S x t ; 20S0 ) , SVS iZi^ , DCK-I L n , i)Kl.M A A 
iwTKnKH tKC2oi>a) , ini 
'XpFiKUflXTS’il ,i)K/ICK='DSk'X 
npF(JC‘J0IT=')i,0CV tCH=:'i)St\’) 

1)11 lf> lal, I 
tlifilaO 

10 conTii'i'Jt: 

lal 

SVCUaO, 1 

RKAi>Ct>l,^)SA1Pbti 

IFC;>AMPbF.f;T.SVCU3(;0 Tu lOX 

CfJ fO 200 
100 IRl»l 

200 CUHTtNOK 



1 l7 


3U0 


400 

500 


600 

7 U 0 


2 U 


10 


C' 

c 

c 

c 

c 

c 

c 

c 


A - 1 . 6 

U — 0 , 6 j > (» o 

r,^(i 3 =;>vci) + giji'uu 
i = r + i 

Tr' (oi'i.i; t’,o«;u -i/tO A ^suv'IjOax 

HEADCbl ),JAMPuh 

TF(bA.v,p,,F,,.'llT,£i\r( l’-U>Gn TO 40U 

IHCI )=•! 

no I'n 6(>f> 

tiUDsI 

COfil'tuJi. 

tFCxotT. ro 6v«) 

nj=-‘urH> 
r,i) T'l 7 «hj 

rut~£«ui) 

r./C n=ovf i-i) + io 

Hi'.' Ci1 

rn i'CTi!'-!; 

rt,, .)( 1-1,1 

:-iJ( L)=o\!U)*^rn'c^ 

••!.<TTP:(6i, + )svtn 
ruiiiHunF 
OFTU'Si 
F. |0 


FtirOintjTcHF apcmco:- 

THTiJ Ct'iK PFFFOhM.S A"CF F-jCoOH'JU i)'' fn’< n.iTP'U 

OF IMIF UECl.lATnu, 
tilPUT;SliB-OAf)u SI'JOAIj. 

n.jT'HU-rscooKu u iu-oajO) oiouai,. 


SIIOKlnJt lOiC APCMCOtOp) 
ruianfj/rt7/BiT,ni:h*nfi 

HtiAL HnoTC4,U;8} ,UKLMlM,l)kiTjMAX,VAO U»I 

tOTEGER YC4) ,Ort,.iP,rTOO 

nilf iT ( 2 , 0 > =0 . « 5 yM? J uT c 2 , t 3 = I , 0 

Mijr/2(J,u)= >,H455MUETCl,l J = l,0?Mnu'rt3,23~l,0;>,(;urt 3, 3)si,n 
M Jl/f C 4 , 0 ) SO . 9 ; MUu f C 4 , I ) =0 . 0 ? HUET (4,23=0.0 ; MOl/l (4,3) =0 , 9 
M.JEX (4,43 = 1.2; MOOT (4,53 = 1 , 6; MuET (4,33 = 2.) ; MOT, i ( * , H = 2 . 4 

3IT3CWt,(wP) 

!>FEMXri = tfAL4IN(BP) 

nEbfiAXsuFuniM^^ioo 

Si s8 0 trf I j M X 13 I 

OPF4(UJ3IT»61 ,0EVtCK:=*!)SK') 

maoc6i,i'3'J 

QPBO ( U.i IT* 62 , DEV ICFs ' DSK ' J 

MRirE( 62 ,*)N 

DU 10 1 = 1, !< 

IFCb.oT.DEEHIN) S=DEIiHIw 



n rj n o n n n o 


1 i 


tK C o , G r , OGL - 1 A 4 ■) Hu; I mM A K 

'-^GAOCdI Gt-; 

nil jsUruT 

yc.n=u 

IFCoAMnUK.Gt.HAiiGk;) GO Tn 20 
S A ^ P G i:; 3 6 A K p b E - H A ‘ > G !■: 

2Ci n,\nG£=n,viQZ/2 

rfRTTE(62,9^^^) tYtJj , J=1 ,nIX) 
oyp FljnMATClXrAilJ 

'iiirpux=o 
ou 4u ,j=2,Hir 
n:,iTPHi=:2*U‘iTPijT+ytj) 
tU C'.v]'ll«nr; 

Ti^ c V ( 1 ) )un‘i'pj'r=c2+4-(aT r-i) j-in r'\)T-i 

G= (Hi TfUnTPiTj 

10 c'n'raii'U. 

CU lolU'luCl’Sol ,r)hVlf!!;=»i).iK'') 

ri.OaE £ ii. r , OoV tCi,:= * u;>k » i 

* 

Gjn 

Gl/HiGIUTlUE APC'IDECHP);- 

TU£.> SUAHOUrilE PERP'JRhR ARC i iJFC’Uit,>*G , 

I,4PiiT:CtjDE0 SUB-BA'llj SIG’MAG. 

a.rrpuT:(»Ecai,)Eu 3na-;iA<4i) sig.iau, 

S li A K J U T 1 ,>1 fc. AP( j M I,' E t A P ) 

Cf« lMlhi/A6/3TnKV 
CG -mOu/AY/iU-tfOEG !l!l 
,''u'^A3A//vB/''»!J ,AYAu 
1)1 -lEUGIUM AVAi*(‘>,2'Jt>',0 

RUAb finurci,o;p) r i/AbUE,HF.AO,VAn,r’ tnEy,nLbNi>j,i)i' i/' \x 
TllTEGER Y£1),BIT,0P 
MllLl'(2,U)s:U,'-15f»rJUTC?,ll = l.‘i 

Aura'(3,U)sO,AA5bMnijT(3,n = l,o;M!JuTt3,2 j = l .u? n>!liT( 3 , 1 ) = 1 . i 
Milbl C i,0)s0.9?MUT/P(4,i)sO,Q;MUbTCA,2)=').9?Mul.i.’('i,3l=0,<^ 
AuTiT( 4^4) = 1 .2; 3UT/fC4,r>) = l-fi?ntib'CC't,H)=2.0 7Pub C f i , 7 t =2 , A 
nKbMAX=i.)Eli3t‘Ktf,-'J 
ssoi'X.irn 

npEo ( a:J iT»B2 , us y ICE= 'dsk * j 

npEfiCi|fUt=H3,UFVICE='USl<'l 
nSADCo2,*)NlW 
i4RXTE(63#*)U» 

no 10 K=l,MW 

1 F ( S , lit , UEbH TU ) SsUETiM Hi 
TP(S,GT,nSLiUX) SsDELMAX 
REAU(B2,‘?99)CnJ),J=lfBIT5 
999 FUH«aTClX,4Xl) 

IVAbUEsO 



f 


0 .) r= 2 ,.ur 

rU(vi 4 >liJ=C'/hij‘j>r 1 + 2 + Y(.U 
20 Cllilini. 

ft'f Y c I ).r:.l, 0 ) lVAt-U^l=fi»+(jiTl-l) )-jVaI/!K '-1 

V40'iK='Hi;VAh'JK-,-:>/2 
IK i ( I ) .n-;.o) 

VAraii:=:v ibOR:*3’cn,-.v 
WKX'i'Ktfi.i , »•) YAIjUl 
AVAb(‘iPKJ=VAOiir. 

1 '=() 

Dn jn IS2/-5II' 

noTt''0'i:=2i‘i)irrn'<jt+Yt>T) 

30 Ciju "u ii: 

CKCYCi ) . lli:.l ).J0n’uTi:C2'»‘*faTT^n Jx.i U'P.iT-t 
As-, ct>ii D 

IJJ niill, M'li-: 

f;j,o..!:u^in'=n2,n.*:viC!:='oAK *) 

Cu )o'K Ji'iET^o i 

it) 

c”"”" 

C rriirtCC):- 

c r.ir.s f'l.aci'iA'i PiinvtuRb 'ffU:. ni^ riir?; :i.jr:)3 

C 'iuC-jniUb to raNuk. 

c 


iiiTi:ai:’-{ Ci’Da,bP 

IKC*)P,K 9 ,t 1 CTOR =3 

TKCiio.:-: 2 . 2 )CTnB =2 

IFfOP.!':-,),4)CTilB=2 
Tf'(0P.K,),S)CTuR = 2 

OtIT.IKtl 


c. 

c 

c 

c 

c 

c 


F.-iO 


'uflCTin i VAuMXNCBP) 

Tills fUaCTIU^J PHnVIDES THli YAbUFS Of UFbHiM Fyn 'KinnjS 
S)JR-B,\.il)S, 


FUHCTtt)'^ 

laTKajiR DP 

Tf C BP , Eg . 1) V AL)M 1 0=0 ,01125 
If CUP,Kg.2)VALMIU=0,Oi3 
If CoP,R0,i)VAl4MX0=U,015 

If CaP,Fg,4)VALMIN=O.Ol5S 
IF ( tiP . HiJ . 5 ) tf 10=0 ,0 1 8 

RETlIRfl 
Fiin 



r 

c t-’.nc’L'i j STsrvj';:" 

c riKJCfiu.j i>Pt)vinn.s T;,T'nAi. 

c 

nuAii 

inriXtiw A(> 

I ?’ ( a p , n y , 1) ST H T «= (' , u i ••> i I) 
tKfaP.!: j.<!)STsc^K=it,na2 

Tr’ C HP , »’.) , i ) STS CZ'':=:a . 002S 

T K ( a p . I / i . :■) ■) a r.'i I S''=a . P') j -j 

■It.Ta'Ti 

r-na) 


c. 

<: 

c 

c 

c: 

c 

c 


r. iasi'aM ii; SPA r.a (B.>) 

i:us sjAKOuTi 'Ji. [fi (!iinoucTif<.. •,! I’M suBa-i.iT t,,'!-. 
U'TKM.-’.i'.ATho T,il' SpF.hCH SA iPjiTo YO riiii: PHlAl iAt4 

a:» PaTr-;, 


SUBKOtJ I’lNt: SRATkIURH) 

C 0 Mfinii/A 2 /rutiaT 
CO'1Tr)H/A5/MXirPC3L,P, taAM») 
i/AU/fltYT 

r>i n;isTij'i xc?A,5h) 

laf’ir.FP !!MAfiKb,aP, COUNT, MX, H(Ujn,Tr>Aan,!U'it»FT.yl 
treat. IT. h:^,o)Gn co i.ju 
rC!ai 521 
C.u TH lit 
100 PUi‘t=X< 

110 COUTCilUlC 

DPHfJ ( u'l ITSPI U 4 , UH: V ICl>:= ' URK ' ) 

TF(aPj:j.UGt) TG 150 

nilUUTa^ 

GU TO '2U0 
ISO CLlONTs? 

200 CDHlMlUt, 

itiA 0 f)=np 
NUliH-7 
RE UNO Ft.i'i 

100 PMAO tFIT.4,*) CHANES 

IK tCHANEiu.EQ.BP) GO W 400 
read (FlLlr^lOJ 
GO TO 100 

400 READ (FID4,*) M,D 

CUUNTaCOUNT-1 



1 . 1 


\0 


c 

c 

c 

c 

c 

c 

c- 


20 

10 


.10 


€• 

r 


r„r.u 

(»nT.4, «-,K ji-= i»i CH.i'lii,!, 
Cl-' ,'n ion 

'1 

) 




SiM^u’.iTuTi uH; ;- 

i'iUS SiJ5ii iwri'irJ uftfn.jp mk i '*!-; 

TfHMIT: Siin-t’.AMu .;ax;>'AliS ’]>'■’ Ai.T. Tii?’ T. 

OUTPd CM'J fPiJT Hir I’dK .JUd-iiAiiD CulK-d' . 


j,., o i , TKlI/iul 
nif *.t,i i 'ipt (fi, 

K,1 i 

!>u r = t ,Ki 

yin~u,p 

rn! /it: Is 1,1 

xtl j=4( L) iHint-nU 

CilMiT ii'i, 

rii (n.v,'"'. 

nt'Ci.C U''iTsTFTiril),uKVT(.:!r;='i)‘lK' j 

nil in 

rii'rrhj'i!; 

ri,fii.>(U”.iTT=.iFi(<:n3,nKViCh:='DtiKM 

ncT.iPti 

Fvm 


C SOnKlUiTtfJti Sr-Jh':- 

r Tiur, .sdHHnDTrdii: cAhCui.AThr> tuf snuAu th nmaF nvfF), 

C IflPnF; IdP'JT AhD nUTFlI'l UK ThE SVSl'E.'i, 

C nuTPilTlalGUAL TU HA-ltii IjK ' TrtF Ci)ni;;P T .1 l/P. 

C 


SdftKnuTlNK h'iH 
CUMfin.'i/n'l/SNRdH 

X(2.'Su) ,^(20*10) 
REAIi AUijMER, AUEn'ia,{-i!>lt{IiU 
Addf-lFfts'I 
ADKdMH=n 

0PF:tiCti"IXT=60,l)EVICR='lASK' J 
npKR c i w rr=6 s.oev tCR= ' dsk ' ) 

REAiJ(hO,#)ttl 
READ Chi,*; 1 

DU 10 talrltS 
READ(t>0,*JXCU 





« 


10 


iLtlf Xl.U *♦2 

A t » f >: i < 1 i -I = u i » f- j f- . K f ( Y 1 1) - X Cl 1 ) * * 

rnr<i rriH|.; 

r..l'UJfi=Ajin il u ( AMil'llCH/AuKwMH ) + 1 0 
f'hf laK ( 'Ij4 TT=60 , nir/XCt=: *IUSK ' ) 
ri,<u,.Kcn.ii r=to3,DhvirK='i).sK ') 
pirPlif.'ii 


SiiBKnijTxfjc 

TffiK SuBR'HPfl 'JK (;AuCM1.4Tk;.< 5 .'t’.U. HuR UK .'iU,. 

CUl^KK. ri<r, UK YUa .SF-X^f'/jUl IB IJh UaUT t,A. 

TtiUti P: 1 iP.J’' A.iU UufpUY iiPEKCH F.A”iUM*,:l UP YHr/ (.'UlU':!-’ 


r\ 


nuYUlfT: oiUilAij I’fi r.'uiUK paYi .1 Ut- 

fi.U-UXf UP (Ml-. P.Im; li;:iT=:l2ii BA'iU.i.-'O 


(til, CUiil'!-' 


•^llUrKM'P CUK .iXCBUR 
PU*l*AUii / / fci lisAt'.*'!' 

!< I Ul.'lUl ■) i X C ) , Y( 2(.'5(! ) 

PUAIi n.*u Ij.R, 4t»!-;,Klk,i-;nK.SKU,Tn 

T 1=1. 

r>uf'.u''k.=i-.p 

jf'U 

900 PUP., M A PC IX, 'tn.PlAUK BPKCXKV TriFKr>HBl,U up o’UU JL ,) ' ) 
tiBU •Uil 

901 PU'?, lATC a, '1-1 BPPCTFTKB T.lRr:.9.i‘'n,n ' ) 

ACC k Pi'*, Yu 

Ui‘PM(il1lT=Hv,MH:/TCK='unK' J 
UKKi- ( tif. 1 ■”=!' 3 , UE V ICC= ‘ DfiK * ) 
lU-iAi.'COt),*)'': 

RCAuft.3, *r; 

DU il> 1=0,15 

K(:Au(t,(-,*J(A(.l),.r=;12H*Ttl ,12n*Tf1 2P) 
(YfU),iJ=t2P*l + t ,J2P *1+1201 
nil /u T=t/.H*,r + } ,1 2B*r+i2« 

AMUMPi',=U,P 

Ar>P'UU(=(i,tt 

AltUUP'k=A£i'l'lKR + X( J)**2 
AuPii'BHs/ini; 'i.'lH+CXf vi)-YcJJ)**2 
20 CU?iiT!<(Ut. 

5i4I7s:AI,UG lu( AUU -iHU/AuRwMH) *lu 
JKCc>MH,GT,Trf)kn Tit ju 
SUREli.i' 

r,u xn 10 

3 0 t> bG=SH KSfciG+S 0 F 

I IslT+t 

40 CU^llriBK 

10 CDUI'IMIK 

snrgpg«smrskg/ii 

Cii0SECUNlT=60,DEVICK='nSK'3 

Cl4Gi>RC'l'<TTao3,nEViCp:='USK*) 



o n o ^ n 


1 j ^ 





,s.]Bitn«!Tr,.)K Biua j-u- 

'i’illB r>IIH«;1.iT [,•;(, HKriTS TUK nV''’KAiif. <?•' SIG.iAu 

T'l llui.in !(4TlfJ OF fHF SOH-rtAftii’i 


GunicnuTiOh; prig-jr 
rri<“!riniJ/A9/RWR|)B 
COOOU,j/alG/o>''.KoFO 
RKAb r.'-lKUn^.SX'HGbG 
WtaTFUG,0 F$J 

995 FlIiO'.A'l'C 15 a, 'S.JR RKl’FUH'HJiCi: ' / 

1 11, X, '=^5; =====:^ = =: = = -=*//) 

ah'J'.’I'HB, •’ J 

oyr, riiBG-i’ct ja, ' nvMO'lub sun =*,fir.4,* obv 

1 n X, 'Gna.lKil PAl* SuR do'// 

r-hTi,R'; 

KO!' 



